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interference generated by power lines are gap-type noise 
caused by electric discharges across line hardware and 
corona noise caused by -the par-tial breakdown o-P -the air due 
•to -the high electric Fields around transmission line 
conductors. Using original data, this research has devel- 
oped a parameter based model of gap-type and corona noise 
that allows the Fundamental noise mechanisms to be 
mathematically or physically simulated. An expression For 
the power spectral density (PSD) oF gap-type noise and 
corona is derived. The energy detection problem is 
Formulated, and using analytical results based on the Hall 
model For radio noise, a robust energy detection receiver is 
developed. Tests oF this receiver using actua 1 and 
simulated data are described. 
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I- INTRODUCTION 



A. PURPOSE 

Observations oF radio interference at high -Frequency 
(HF) receiver sites have indicated that existing models o-F 
man-made radio noise are inadequate to describe the observed 
time-, and frequency-domain behavior o-F the noise. The data 
also shows that quite o-Ften a speci-Fic noise source pre- 
dominates in a particular location. These observations led 
to the objectives o-F this dissertation: 

1. to develop speci-Fic models For certain man-made 
noise sources in the HF radio band and 

2. to apply this model to the analysis oF energy 
detection rece i vers . 

B . BACKGROUND 

Studies conducted over the past Few years have indicated 
that one oF the primary sources oF man-made radio noise in 
the HF band are alternating current (AC) transmission and 
distribution lines CReFs. 1,2D. Two oF the primary sources 
oF power line noise are gap noise, also known as micro- 
sparking, and corona. Gap noise is caused by a sparking 

process between hardware points on utility poles and corona 
is caused by the partial breakdown oF air due to high 
electric Fields around high voltage conductors. Both oF 
these noise types are non-Gaussian noise processes. 
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The majority of the models of man-made and atmospheric 
noise to date have assumed the high amplitude impulses 
driving the receiver have either a Poisson arrival rate or a 
variation of Poisson arrivals. This assumption has allowed 
the derivation o-f -first order envelope statistics: the 
amplitude probability distribution (A PD) -function and the 
density -function o-f the phase o-f received noise -for various 
amplitude distributions o-f the driving impulses. Observa- 
tions o-f man-made noise -from power lines, however, has 
indicated that the assumption o-f interpulse independence 
(Poisson arrival times) is not valid. There exists a 
de-Pinite time domain correlation o-f the impulses driving the 
receiver. This is due to the underlying deterministic 
mechanism o-f the -fundamental -frequency o-f the power line 
voltage. Consequently, this research has concentrated on 
the statistics o-f the impulse arrivals which are manifest in 
the autocorrelation -function and the spectrum o-f the 
observed noise process. Accordingly, the models developed 
are speci-Pic rather than general and are more suited to 
source identification by spectrum analysis and robustness 
evaluation o-f systems rather than to generic specification 
of optimum receiver structures. 

The performance of energy detection receivers for 
stochastic signals in non-Gaussian noise has received 
relatively little study. Actual implementation of algo- 
rithms for signal detection has been based on heuristics. 
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The advent ot digital signal processing techniques that 
allow tor complex post-detection algorithms suggests a 
careful statistical analysis of the noise performance of 
these systems may help develop more etticient and robust 
rece i vers . 



C. CONTENT 

A briet outline of the rest of the thesis will now be 
presented. Chapter II develops a general model tor bandpass 
impulsive phenomena using a tiltered impulse model and 
complex envelope theory. Well known empirical models and 
physical models ot radio noise are presented and their 
relationship to the generalized model discussed. Particular 
attention is paid to the Hall model tor atmospheric radio 
noise CRet. 3D. Some usetul extensions to it are derived 
that will be used tor simulations in Chapter V. 

In Chapter III, tield observations trom sources ot gap 
noise interterence are presented and analyzed. Three cases 
ot actual interterence are used and a probabilistic model 
based on a statistical analysis ot the data and the ti ltered 
impulse model developed in Chapter II is specitied. The 
power spectral density (PSD) ot gap type noise as predicted 
by the model is derived and compared to the data. 

Using an actual case ot corona noise trom a 500 kV power 
line, the same type ot analysis is carried out tor corona 
noise. The corona noise model is based on the same tiltered 
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impulse -Framework as -the gap noise model o-P Chapter III 
however, -the specifications di-P-fer substantially. 

Using the Hall model -Por atmospheric noise discussed in 
Chapter II, a locally optimum receiver -Por detection o-P 
unknown signals in HF atmospheric noise is derived. A 
practical mod i -Picat ion to this receiver shows it to be an 
adaptive limiter. The performance o-P this receiver is then 
examined in simulated man-made noise and in recorded HF 
signal and noise data. 

Appendix A presents a brie-P description of the 
instrumentation that was used to collect much o-P the data 
-Por this thesis. The dissertation concludes with a brie-P 



summary o-P results and some suggestions -Por -Future research. 



I I . IMPULSIVE NOISE AND RECEIVER MODEL 



A. INTRODUCTION 

In this chapter a general set of specifications -for an 
impulsive noise interference process will be developed and 
the results o-P previous work in the -field explained in terms 
o-P the specifications. This generalized scenario will be 
used in later chapters to speci-Py complicated sources o-P 
man-made radio interference. In attempting to describe 
atmospheric radio noise, two general types o-P models have 
been developed: empirical models designed to -Pit -first order 
statistical data, and physical models directly related to 
the underlying physical mechanisms. The Hall empirical 
model -for atmospheric radio noise has been shown to -Pit 
atmospheric noise data very well and will be used in this 
dissertation. Some extensions to the Hall model will be 
developed in this chapter. 



B. GENERALIZED SCENARIO 

In order to provide a -framework -for the discussion o-P 
the impulsive noise models, a general interference scenario 
-for impulsive noise will be described CRe-P. 4,53. A typical 
interference scenario consists o-P the following elements: 

1. a source of i nterf erence , 

2. a transmission medium to the receiver and 
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3. the receiver where the interference manifests 
itselF. 

For this study the noise process at the receiver is 
modeled as the sum of a high density (in time), low amp- 
litude Gaussian component and a low density, high amplitude 
impulsive term. 

These .elements are shown below, where e 7 (t) in Fig. 1 is 
the impulsive interference, L 7 (a>) is the Frequency response 
oF the transmission medium, z'(t) is white Gaussian noise 
and H '((tf) is the combined response oF the RF and IF Filters 
oF the receiver. 



z 7 (t) * 



e 




Figure 1. InterFerence Scenario 
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The source oF impulsive noise is specified by the 
equat i on 

N(t) 

e 7 (t) = E e t 6 (t-t 4 ) (2-1 ) 

i = 1 

where e 4 is the real amplitude oF the ith pulse and N(t) is 
a. unit counting process that generates arrival times, the 
ti 7 s. Using this representat i on , the impulsive interference 
source is described by the probability density function 
(PDF) of each e t , given by p t (e) , and the impulse arrival 
times generated by the unit counting process, N(t). 

In order to simpliFy later analysis the interference 
scenario will be expressed in terms of its complex envelope. 
The interference process at the receiver detector is n 7 (t) 
and will be replaced by its complex envelope equivalent, 
n ( t ) , where 

j<Do t 

n 7 ( t ) = Re C n ( t ) e 1 (2-2) 

and 

n ( t ) = n c (t) + jn, (t). (2-3) 

J t 

The term n(t)e is known as the analytic signal 

representat ion oF n 7 (t), and n(t) is the complex envelope 
representat i on o-P n 7 (t). The reference Frequency For the 
complex envelope representat i on is u) 0 . When n 7 (t) is a 

bandpass process, n c (t) and n* (t) are the lowpass inphase 
and quadrature modulation components oF n(t) respect- 
ively. CReF. 6:p. 75D 
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SpeciFying a -Pi ltered impulse model in complex envelope 
form requires deriving an expression for the driving impulse 
-function given in Eqn. 2-1. Two equivalent approaches may 
be taken and they are both outlined below. In the -first 



approach 


a 


bandpass impulse 


may be postulated 


where the 


spectrum 


o-P 


the impulse is 


assumed -Plat From 0) o 


- W to U)o + 


W where 


(do 


is the reference 


Frequency and W is 


an arbitrary 



bandwidth. For this case the analytic signal representat 1 on 
is easily -Pound in the -frequency domain. Using the 
■frequency domain deFinition, the analytic signal associated 
with e'Ct) in Eqn. 2-1 is 



e A (t) 



C 1/n ) 



jut 

E ' (u)) e da) . 



0 



(2-4) 



Inserting the de-Pinition o-P the bandpass impulse de-Pined 
above and using Eqn. 2-1 



N(t) jtiJo (t-tt ) 

e A (t) = E (2/n ) e si n ( W (t-t 4 ) ) / (t-t» ) . (2-5) 

i — 1 

The complex envelope o-P e'(t) is de-Pined as 

- j(do t 

e(t)=e A (t)e (2-6) 



and 



N ( t ) 

e ( t ) = E (2/n 
i = 1 

where u) 0 tt = 0 4 . As 
tically approach CRe-P. 



” jQ i 

e sin(W(t-t t ) ) / (t-t 4 ) (2-7) 



W becomes large, e(t) will asympto- 
73 
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N(t) -J0! 

e(t) = E 2e & (t-t 4 ) . (2-8) 

i = l 

The second approach CRe-F. 83 is to consider "the impulse 
as ideal and constant in the -Frequency domain. The analytic 
signal representation o-F e'(t) is then directly determined 
•From the time domain de-Finition o-F an analytic signal. The 
analytic signal is de-Fined as 

e* C-t) = e'(-t) + je'(-t) (2-9) 

where e'(t) is the Hilbert transform o-F e'(t). The analytic 
signal representat i on o-F Eqn. 2-1 using the above defin- 
itions is 

N(t) 

e* (t) = E 6 (t-tt ) + j/(n(t-ti)) (2-10) 

i = 1 

since 1/nt is the Hilbert transform o-F the delta -Function. 
The complex envelope o-F e'(t) is then 

N(t) - J0 t -JCDot 

e ( t ) = E <S (t-ti ) e + je /(n(t-t 4 )) (2-11) 

i = l 

using Eqn. 2-6. 

Although apparently di-F-Ferent, when the complex impulse 
trains described by Eqns. 2-8 and 2-11 are convolved with a 
complex -Filter impulse response, the -Filter output is 
identical -For both representat i ons . I-F h(t) is the complex 
envelope o-F the -Fi Iter impulse response, then the f l Iter 
output with the input given by Eqn. 2-8 is 
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-J01 



( 2 - 12 ) 



N(t) 

E e h(t-t 4 ) 
i = l 

where complex convolution is defined as 



h ( t ) *■ e ( t ) 



( 1 / 2 ) 



h(T)e(t-T)d , T 



( 2 - 13 ) 



The f i I ter output with the input defined by Eqn. 2-11 is 

N ( t ) - j0 i - jo) 0 t 

(1/2) E h ( t-ti ) e + je h A (t-ti ) (2-14) 

i = l 

where the convolution definition of the Hilbert transform 
was used CReF. 6:p. 69]. Using Eqn. 2-6, the same result as 
derived in Eqn. 2-12 is then -Pound. Due to its simplicity 
and for ease o-P use in programming simulations, Eqn. 2-8 
will be used to represent the driving impulse function for 
•the rest o-P this dissertation. 

One further point to mention concerns the probability 
distribution o-P 0 t in Eqn. 2-8. For the processes we con- 
sider, t 4 will have random arrival times and, as can be seen 
in Eqn. 2-7, 0 t is formed by multiplying t t times the 
reference Frequency, o) 0 . Consequently, 0 4 is the phase of 
the impulse time with respect to the reference Frequency. 
It can be assumed to be uniFormly distributed over 0 to 2n 
when the t t 's have a probab l 1 l st l c interarrival distribution 
and u)o is much greater than the inverse oF the interarrival 
times CReF. 9: pp . 279, 10]. 
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In an actual interference situation -the received 
impulses are f 1 ltered by the radiating antenna, spreading 
losses, atmospheric attenuation, receiving antenna and 
cabling losses. These -berms will be lumped together as an 
equivalent -Filter, L'Cw). Examples o-F this -Filtering term 
-For various interference sources have been recorded by many 
researchers CReF. 11,123- The -Fi Iter char acter 1 st 1 c is 
specified in the -Frequency domain as L'Cuj) since that is 
where it is most easily observed. 

An example o-F this type o-F -Filtering is shown in Fig. 2 
where the -Frequency spectrum o-F an impulsive noise source 
due to a gap discharge is shown -From 0 to 200 MHz. The 
straight lines in the 3-Axis view are due to stations in the 
high Frequency (HF) band From 2 to 30 MHz and the FM band 
From 88 to 106 MHz. The continuous envelope seen in the 

upper view is the magnitude oF the Frequency response oF 
L'Cqj). The solid line in the upper picture is the noise 
Floor oF the spectrum analyzer at -100 dBm. 

The complex envelope oF the impulse response oF the 
attenuation term, L'(u)), is given in the Frequency domain by 
L(u)) = Lp CL'C-ai-too ) + L'" (u)+Q) 0 3 (2-15) 

where Lp denotes the lowpass part oF the quantity in 

brackets . 

A complex Gaussian component z(t) will also exist in the 
observed noise process and is due to the combination oF: 
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Figure 2. 



Gap Discharge in Frequency Domain 



dBm 
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1. "the thermal noise in the receiver, and 



2. the combined sum of many low level atmospheric or 
man-made impulsive sources. 

This term will be modeled by independent zero mean Gaussian 
noise processes -For the inphase and quadrature terms with a 
power spectral density (PSD) in each component o-F N 0 /2. The 
distinction between the Gaussian noise term and impulsive 
noise term is that many receiver responses overlap tor the 
high density, low amplitude Gaussian case, satisfying the 
condition for the Central Limit Theorem. For the impulsive 
noise the receiver response to each impulse is discernible 
such that the probability of more than one or two pulses 
overlapping is negligible. Common sources of impulsive 
noise are man-made noise due to power lines and atmospheric 
noise due to lightning. 

The time domain input to the receiver is now modeled by 
the expression 



n ( t ) 



(h(t)/2) * [ z ( t ) 



N ( t ) - j 8 

+ E 2e t e 
i = l 



CD 

j cu ( t-t t ) 

e L ( w ) d(j 1 



(2-16) 



where 

h ( t ) = h c (t ) + j h 3 ( t ) . 

The noise power bandwidth of H ( o> ) will be designated B e F f . 
In most situations L'(oj) will effectively be constant in 
comparison to the narrowband receiver filter H(u). This 
makes L ( oj ) a constant that depends only on w 0 . 

Add i t i onal ly , g(t) will be defined by 
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g (t ) 



(h(t/2) ) * z (t) 



(2-17) 



and is "the low level Gaussian noise at the detector. The 
power in the inphase and quadrature components oF the 
Gaussian noise will be N©B«»tF/2. The complex lowpass noise 
process can now be written as 



N(t) -j0 t 

n ( t ) = g ( t ) + E a, e h(t-t t ) (2-18) 

i = l 



where 



a* = ei L ' ( cdo ) . (2-19) 
The parameter a« is then a scenario dependent parameter that 
is a function of* the reference frequency and the impulse 
amplitude. With the above simplifications, Fig. 1 can now 
be expressed in complex envelope Form as shown in Fig. 3. 



z (t) 




Figure 3. Complex Envelope of Interference Scenario 



The terms From Fig. 1 have been replaced by their 
complex envelope representat ions and the prime is 
indicate this. 



equ i va 1 ent 
dropped to 
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Some further simplifications will now be considered 
based on the character istics oF the H'(cd), -the bandpass 
fi Iter. I F the impulse response h(t) oF the -Filter is real, 
H'((d) has conjugate symmetry about the origin in the sense 
that H ((d) = H ' ( - id) . Furthermore , i F H ' (to) is sym-metnca 1 

in the same way about the reference Frequency, (Do , then 
h s (t) will equal 0. 

The envelope squared of the process defined in Eqn. 2-18 
l s 

2 2 2 

E (t) = n c (t) + n 3 (t). (2-20) 

IF the Gaussian noise term is assumed to equal 0 then 
2 N ( T ) N ( T ) 

E (t) = E E a t a 4 h c (t-t* ) h c (t-tj )cos(0 t -0j ) + 
i=l j=l 

N ( T ) N(T) 

E E a t aj h 9 (t-t t )h 8 (t-tj )cos(0i -0j ) . (2-21) 

i=l J=1 

IF N(T) deFines a low density counting process such that 
there is negligible overlap between subsequent pulses and 
h(t) is the complex envelope oF the impulse response oF a 
bandpass symmetrical Filter, then the envelope oF n(t) can 
be simply written as 

N(t ) 

E(t) = E at h(t-ti ) (2-22) 

i = l 

where h(t) is equal to h c (t) . 

At this point oF the analysis two avenues may be 
pursued. IF there is a time dependency between pulses (non- 
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exponent: i a 1 1 y distributed interarrival times), 



then tne 



correlation -Function and spectrum o-f the process can be 
examined. It the amplitude probability distribution (APD) 
o-f the envelope is the desired result, then the char- 
acteristic -Function o-P the envelope is most useful . The APD 
-Function is commonly used in radio noise research and is 
equal to one minus the cumulative distribution -Function 
(CDF). In general these two approaches are mutually 
exclusive since the assumption o-P time dependency between 
the pulses makes the envelope APD -Function di-F-Ficult to 
calculate. I -F the impulse arrival times -Form a Poisson 
point process then the envelope APD may be calculated: 
however, the spectrum is constant. The remainder o-f this 
chapter will -Follow the envelope APD approach. Results 
obtained by previous researchers will be developed using the 
general noise model. Chapters III and IV will explain the 
other avenue and look at the spectrum o-P complicated man- 
made noise processes. 

C. EMPIRICAL MODELS 



Empirical models o-f atmospheric noise have been 
developed to provide a mathematical expression -For the 
-First-order statistics o-P the envelope o-P the received 
wave-form. In particular, the APD has been emphasized. This 
type o-F model attempts to construct a mathematical express- 
ion that -Fits observed data without regard -For the physics 
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oF the interference scenario. 



Empirical models have the 



advantage that the resulting expressions are much simpler 
than those obtained -Prom the -Filtered impulse models. One 
disadvantage oF this class o-P model is that since only the 
-First order statistics o-F the noise are considered in 
developing and -Fitting the model, the higher order sta- 
tistics may not match the data well. 

Variations o-F the Rayleigh distribution have been 
proposed by a number o-F authors to -Fit observed atmospheric 
noise data. The Rayleigh PDF is 

2 

-a r 

p(r) = 2are r 

= 0 r 



2 0 (2-23) 

< 0 



and the APD or exceedance probability 

2 

-ar 0 

Pr(r > r a ) = 1 - P(r 0 > = e r 0 > 0 (2-24) 

=1 r 0 < 0 

How well the single parameter Rayleigh distribution Fits 
observed atmospheric noise data can be seen in Fig. 4. The 
data points are exceedance probabilities plotted in decibels 
(dB) above the root-mean-squared (RMS) value of the received 
envelope. The data was measured on an ARN-2 receiver at a 
center Frequency oF 10 MHz at Boulder, Colorado CReF. 131. 
The high probability, low amplitude portions oF the observed 
data approaches the slope oF the Rayleigh distribution curve 
For a equal to 12. (a is deFined in Eqn. 2-24) This value 
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Figure 4. Atmospheric Noise and Rayleigh Distribution 
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o-F a was chosen "to -Fit the high probability portion o-P the 
data. That the data approaches a Rayleigh distribution is 
due to the e-P-Pect o-P the many small overlapping impulses 
occurring at low amplitudes whose quadrature Gaussian 
components will have a Rayleigh envelope. However, at the 
low probability, high amplitude part o-P the curve the 
Rayleigh distribution predicts -Par -Fewer values than 
actually occur. This is the distinguishing -Factor o-P the 
envelope distribution o-P impulsive noise: that its 
distribution contains a higher probability o-P high amplitude 
terms than is predicted by the Rayleigh distribution. The 
second curve in Fig. 4 where a is equal to .5 is the 
Rayleigh distribution with the same power as the observed 
data. It can be seen that this curve under estimates the 
high amplitude and over estimates the low amplitude values 
o-P the observed distribution. It should be noted that this 
type o-P plot emphasizes the low probability portion o-P 
distribution since that is where the deviation -Prom the 
expected behavior occurs. 

In an attempt to correct the poor -Pit o-P the Rayleigh 
distribution at high amplitude levels, Likhter CRe-P. 143 
proposed a combination o-P two Rayleigh distributions: 

2 2 
-a r 0 -Br 0 

Pr ( r > r 0 ) = (1 - c)e + ce . (2-25) 
This -Formula has been shown to agree poorly with practical 
results CRe-P. 153. 
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Spaulding, Roubique and Crichlow CReF. 163 combined “the 
Rayleigh distribution with a “power" Rayleigh distribution 
"to obtain a distribution that -Pit very well For atmospheric 
noise over a wide range o-P receiver bandwidths. This APD is 

2 

-ar 0 



Pr C r 


> r 0 


) = e 


r 0 


< ft 
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1/s 










-(ar 0 ) 








Pr ( r 


> r 0 


) = e 


r 0 


2 ft 


(2-26) 


where a, ft and s 


are 


determi ned 


From 


measured 


stat i s"t i c a 1 


parameters oF the 


no i se 


m 









Horner and Harwood CReF. 17] proposed the two parameter 
log-normal distribution and -Pound it gave a satisfactory Pit 
to r^^dio noise data in the VLF band. The log-normal was 
chosen since it has a more impulsive tail and -Pits the high 
amplitude, low probability data better than the Rayleigh 
di str i but i on . The PDF For the log-normal distribution is 

2 2 

1 - ( 1 n ( r ) - a) /2a 

p(r) — e (2-27) 

1/2 

ra(2n) 

Two examples of the APD of the log-normal distribution are 
shown in Fig. 5 where the curve with a equal to -.24 and a 2 
equal to .24 is the log-normal distribution with the same 
mean and mean square parameters as the curve with a equal to 
.5 in Fig. 4. Comparing the two curves the log-normal 
distribution provides an intrinsically better Fit to 
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Percentage o-f Time r 0 is Exceeded 



Figure 5, Atmospheric Noise and Log-Normal Distribution 
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impulsive noise data at: both the high and low amplitude 



values. When the log-normal density has both parameters 
chosen to match the mean and mean square values o-F the 
observed data an excellent -Fit to this particular data is 
obtained. This is seen in the curve with a equal to -1.4 
and <7* equal to 1.4. 

Another model that has been successfully used to -Fit the 
observed APD o-F atmospheric radio noise is Hall's general- 
ized H t M distribution CRe-F. 33. This model is unique in 
that it is not o-F the -Fi ltered impulse type; however, in 
contrast to the strictly empirical models it does postulate 
a random process. 

One o-F the problems o-F Gaussian models o-F atmospheric 
noise is that the dynamic range o-F the model is less than 
the observed dynamic range o-F the noise. To achieve a 
greater dynamic range starting with a Gaussian process, Hall 
proposed a model which considered the received noise to be 
o-F the -Form 

n ( t ) = a (t ) s (t ) (2-28) 



where a(t) is a slowly varying stationary random process and 
s(t) is an independent, narrowband Gaussian process. An 
analytically tractable distribution was selected -For a(t) 
which was chosen to give good agreement between the model 
and measured atmospheric noise data. The distribution -For 
a (t ) is: 
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(m/2) 



1 



2 2 

- ( m/2a <r ) 



p( a) = 

m 



— — - e 
m+1 



<7 r (m/2) I a | (2-29) 
where m and <r are -the two parameters used to tit the model 
to the data. The distribution of s(t) is Gaussian; 



2 2 

1 -(s /2cr t ) 

p ( s ) = e 

2 1/2 



Hall calculated 



p ( n) 



where 



(2n<r i ) 

the density -Function oF n(t) to be 

0-1 

r (9 /2) y i 



1/2 2 2 0/2 
r ( (0-1 ) /2 ) n (n + y ) 

1/2 

y = m (at /a) , 



0 = m = 1 , 



(2-30) 



(2-31 ) 



and T is the gamma -Function. For the case where at = <s , 
Eqn. 2-31 is the density Function For Student's "t" 
distribution which is the basis For describing the density 
as a generalized "t" distribution. 

Hall then calculated the envelope distribution For n(t) 
based on the above assumptions and obtained 

0-1 r 

p ( r ) = (0 - 1) y 

2 2 ((0+l)/2) 

(r + y > (2-32) 

with the phase uniFormly distributed between 0 and 2n . 
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For atmospheric radio noise, values of 0 in “the range 2 
-to 5 have been Found -to give excellent agreement with the 
data CReF. 18,193. (Models with integer values oF the 
parameter 0 will be abbreviated to Hal 12, Hal 13, etc.) 

Fig. 6 shows the Hal 13 Fit to the same atmospheric noise 
data used in Figs. 3 and 4. An excellent Fit to this data 
is obtained. The Hall model has some disadvantages. The 
higher moments only exist For orders greater than 0-1 and 
the parameters 0 and y must be determined For each 
interFerence scenario and are not easily related to the 
physical source oF the interFerence. 

D. FILTERED IMPULSE NOISE MODELS 

The Filtered impulse models diFFer Fundamenta 1 ly in 
concept From the empirical models described previously. 
This class oF model will be explained in terms oF the 
generalized interFerence scenario presented at the beginning 
oF the chapter. The Fi ltered impulse models have the 
advantage oF being based on the underlying physics oF the 
process but suFFer From the disadvantage oF being 
analytically complex. Interestingly, some oF the results 
obtained From the Filtered impulse models have retroact i ve ly 
justiFied the expressions that were derived From the 
empirical models. 

To determine the APD oF the envelope oF a Fi ltered 
impulse process, the joint character i st i c Function oF the 
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Percentage o-F Time r 0 is Exceeded 



Figure 6. Atmospheric Noise and Hall Distribution 
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inphase and quadrature components is determined. Since the 
noise is narrowband it will be circularly symmetrical and 
the 2-D Fourier transform will become a Hankel transform 
CRef. 20,211. This character i st 1 c function is transformed 
to polar coordinates, then inverted to determine the APD . 
Starting with the complex noise process described in Eqn. 2- 
18, the joint character i st i c function of the inphase and 
quadrature components is 



j (u)i n c (t) + u) a n s Ct) ) 

0 (oil ,(i)a ) = ECe ] . (2-33) 

It is shown in Appendix B that the char act er i st l c function 
of the envelope due to the Gaussian term is 



-NoBCi d„ )/2 

0a (u>r ) = e (2-34) 

and the character i st i c function of the envelope due to the 
impulsive term is 



X 



p ( a ) 



C Jo (a), ah c (t)) - 1 ]dTda 



(2-35) 



0 



0 



0 1 ( (D r ) = e 



where p(a) is the impulse amplitude density function, X, the 
Poisson rate function, h c (t), the inphase component of the 
filter impulse response and J 0 is the ordinary Bessel 
function of the first kind. The density function of the 
envelope is then the inverse Hankel transform given by 



p ( r ) 



CD 

r(i) f Jo (r U) r )0<a ( (Dr )0I ((D r ) d(Dr . 



(2-36) 
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This equation for the PDF is diFFicult to evaluate; however, 
a number oF researchers have examined special cases to 
obtain results. 

Furutsu and Ishida CReF. 223 derived an equation of the 
same Form as Eqn. 2-36 with no Gaussian term. They obtained 
many approximate results by considering the equation For 
speciFic source density Functions. IF p(a) is exponential 
they Found that the resulting envelope density was 
approximately Rayleigh For small amplitude values. IF p(a) 
was uniFormly distributed, the log-normal density was Found 
to be a good approximation over most oF the range oF 
interest. IF a strong local source oF interFerence was 
present, Furutsu and Ishida showed that a Function involving 
the conFluent hypergeometr i c Function was a good 

approximation. They also considered one special case were 
the unit counting Function N(T) was modiFied to be a Poisson 
Poisson branching process and showed that For low 
amplitudes the resulting density Function was approximated 
by the simple Poisson case. 

Giordano CReF. 233 evaluated Eqn. 2-36 For various 
distributions oF p(a) which he determined From assumed 
propagation laws and spatial distributions oF impulse 
sources. In one particular case he adopted the Following 
assumptions : 

1. No Gaussian component 
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2. Uni-Form spatial di str i but i on of impulse sources 

3. Inverse distance law -For received -Field strength 

4. Arbitrary receiver envelope response 

and showed that the APD o-F the envelope -Function was approx- 
imated by 

K 

Pr(r > r 0 ) = — — (2-37) 

2 2 1/2 
(r 0 + K ) 



where 



T 



K 




h (t ) dt 



j 



0 



T is the observation interval, h(t) the -Filter response, and 
c is a scenario dependent constant.- This result is sig- 
nificant because it has the same -Form as the APD o-F the 
Hal 12 model and physically justifies what had heretofore 
been an empirical model. 

Middleton CRef. 24D included the Gaussian noise term in 
Eqn. 2-36 and expanded the character ist ic function without 
taking the expectation. Then by inverting the char- 

acteristic function term by term, a canonical form of the 
envelope density function can be obtained. The exceedance 
probability was then shown to be 
Pr ( r > r d ) = 

m 

2 2 a> ( - 1 ) m moc ma 

-ro Cl - r 0 E A. TCI + ) l F l (1 ; 2; r 0 ) 3 

m=0 m ! 2 2 

e 

(2-38) 
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where 4 F| is a confluent hypergeometric Function, and a and 
A. are two parameters that are determined by the source dis- 
tributions and the propagation law. 



E. HALL MODEL EXTENSIONS 

For later work in this dissertation a model For 
atmospheric noise statistics is required. This section will 
Further describe the Hall model as applied to HF atmospheric 
noise and will derive some extensions to it that will be 
used to digitally simulate atmospheric radio noise. 

The Hall model For the envelope statistics oF radio 



noise is 



0-1 



p(r) = (0 - 1) y 



2 2 ((0-D/2) 

(r + y ) 



(2-39) 

and was chosen For a number oF reasons . The most important 
are: the Hall model is analytically tractable and shows a 
good Fit to the data. It can be obtained For certain values 
oF the parameter 0, starting From the physically based 
Filtered impulse models. This was shown by both Giordano 
CReF. 23] and Middleton CReF. 24] . Additionally, SchonhoFF 
CReF. 18] has related the Hall models to CCIR Report 332 
CReF. 25], the standard reFerence oF First order statistics 
oF atmospheric noise. In the CCIR report the APD curves are 
shown as a Function oF V 0 . The parameter , V 0 , is a commonly 
used measure oF the impulsiveness oF atmospheric radio noise 
and is deFined as 
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2 1/2 
ECr ] 

V D = 20 log 

ECr] (2-40) 

where r is the amplitude of the received envelope. For 

quadrature Gaussian noise V D = 1.05 dB and as the noise 

envelope becomes more impulsive, V D will increase. 

One problem in relating the Hall model to V D is that the 
moments of the Hall model only exist For orders greater than 
0-1. This means that For the Hall2 and Ha 1 1 3 models, V D is 
undeFined unless the model is modiFied. Two methods have 
been proposed to do this. Hall 's method CReF. 3] adds 
another term to the envelope density Function and multiplies 
by a negative exponential term to reduce the tail oF the 
density Function and give Finite moments. The advantage oF 
this method is that the density Function exists From 0 to 
inFinity; however, the resulting density Function is 
complicated. The method proposed by SchonhoFF CReF. 18] 
truncates the Hal 12 and Hal 13 d i str i but l ons above a set 
level, T P , and renormalizes the density Function by picking 
k such that 

T P 

-i 

k p(r)dr = 1. (2-41) 

V 

0 

The level, T P , is chosen to set a desired V D For the 
distribution. The constant k will be a Function oF D where 

2 1/2 

D = (1 + ( T P /y) ) . (2-42) 



37 



The existence o-F the -First and second moments -For these two 



distributions is then assured. Using this procedure 
Schonhoff generated a -Family o-F distributions parameter i zed 
by 0 that can be used to represent a wide range o-F measured 
data. By matching the V D ratio o-F the data to the dis- 
tribution, a close -Fit to the APD is -Found. 

An extension to the basic Hall model to be used in a 
later chapter is the density -Function o-F a constant 
amplitude sinusoid plus the assumed noise. This density is 
analogous to the Rician density in Gaussian noise theory. 
It will be used in simulating the performance o-F energy 
detection receivers -For. the signal plus noise case. It was 
derived in CRe-F. 26] -For the Hal 12 model and it is derived 
-For the Hal 13 and Hal 15 models in Appendix C. For the Hal 13 
model the density -Function is 



2 2 2 2 
tr+A+y ) 2y rD 



p ( r ) = 



4 22 22 22 242 

(r - 2A r + 2A y 2+ 2y r 2+ A + y )(D- 1) 



(2-43) 



where A is the signal amplitude, y is a Hal 13 noise 
parameter and D is de-Fined in Eqn. 2-42. Likewise the 

density -Function -For the Hal 15 model is 



where 



p ( r ) = 



4 2 2 

2y r ( b + 2a ) 

2 2 5/2 

(a - b ) 



(2-44) 



a = r 
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and 



b = -2Ar. 

Fig, 7 is an example of the APD For a sinusoid plus Hall3 
noise with the RMS value of the noise set to 1 and the 
signal to noise ratio at ”26, -6 and 14 dB . For the -26 dB 
SNR case the APD is very similar to the Hal 13 noise only 
density shown in Fig. 6. 

Commonly used Functions relating to the Hall models are 
shown in Table I. These are the envelope density Function, 
the 1st and 2nd moments oF the envelope, the inverse oF the 
CDF, the quadrature marginal density Function and the Hall 
noise plus random sinusoid density. The parameter D in 
Table I is deFined by Eqn. 2-42 and the Function EC* 3 is the 
complete elliptic integral oF the second kind. The quad- 
rature components can be derived by transForming the 
envelope density to rectangular co-ordinates and integrating 
with respect to one component. The inverse oF the envelope 
CDF is also included since it provides an easy method to 
generate random deviates. 

Also useFul For simulation purposes are the values oF 
the Hall parameters , 0 and y, For diFFerent values oF V D 
normalized to a unit root mean square value. The truncation 
point Tp was Found For a des i red V D using an iterative tech- 
nique. These parameters are shown in Table II. 
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Figure 7. Hall Noise plus Signal Distribution 
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TABLE 



Hall No 



EC n (t)U ECn (t) D T P 



2.1 


5 


1 


.785 


— 


1 .01 


— 


3 
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1 


.708 
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.708 
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.631 


11.66 


.420 


27 . 7 & 
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.562 


25.90 


.364 


71.05 


6 


3 


1 


.501 


67.71 


.321 


211.1 


8 


2 


1 


.398 


11.77 


.085 


137.6 


10 


2 


1 


.316 


17.08 


.058 


291 .0 


12 


2 


1 


.251 


24.25 


.041 


587.1 


14 


2 


1 


. 199 


33.81 


.029 


1141.9 



TABLE II 



F . SUMMARY 

Both empirical and analytical Filtered impulse models 
have been successful in explaining the First order envelope 
statistics oF bandpass radio noise processes. However, in 
obtaining these results the assumption has to be made that 
the impulse arrival times due to the noise process Form a 
Poisson point process. In general this will not be true For 
man-made impulsive noise and For some types oF atmospheric 
noise CReF. 27D. In most cases oF man-made noise the 
underlying deterministic mechanism oF the device causing the 
noise will modulate or regulate the impulse distribution. 
The Hall model For atmospheric noise has been examined in 
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detail and some useful extensions to it have been derived 
that will allow noise and signal plus noise environments to 



be simulated 



III. GAP NOISE MODEL 



A. INTRODUCTION 

In this chapter a speciFic source oF man-made radio 
noise will be modeled and speciFied in terms oF the counting 
process, N(t), and the probability density function (PDF) o-P 
the impulses, p(a). These specifications were introduced in 
the previous chapter to de-Pine the impulsive noise model. 
The noise source that will be examined is gap type 
discharges and is commonly -Pound on electric power dis- 
tribution and transmission lines. Gap noise is one o-P the 
major types o-P interference -Prom power lines and is 
Frequently observed as the primary interference to com- 
munication systems operating in the high Frequency (HF) 
radio band. 

B. GAP NOISE MECHANISM 

At least two mechanisms have been Found by which a gap 
discharge process can occur on a power line. The resistance 
in the line insulators can be degraded allowing current to 
Flow through the insulator base and creating a potential 
gradient across any gaps or detects in the insulator 
mounting hardware. A second way in which a potential can be 
created across an air gap is by an electro-static coupling 
oF the line potential to isolated hardware on the pole. In 
both cases the potential across the gap is discharged by the 
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voltage breakdown o-F the gap and the resulting rapid current 
-Flow or spark. This process generates a radio frequency 
(RF) noise impulse with spectral components extending into 
the hundreds o-F MHz. During a single discharge, the 
potential across the gap is temporarily diminished. How- 
ever, while the -Fundamental 60-Hz waveform is still greater 
than a threshold voltage, the process can occur again. The 
spark will discharge across the gap . repeatedly until the 
alternating current waveform drops below the breakdown 
threshold potential. CRe-F. 28 : p . 78,29] 

The spark discharge and recharging o-F the gap potential 
indicates that this type o-F process is regenerative and can 
be modeled as a renewal process where the renewal points are 
associated with the sparks. One -Feature o-F gap noise that 
complicates the modeling is the 120-Hz on-off-on modulation 
imposed on the renewal points by the alternating current 
<AC) waveform. One way to account -For this e-F-Fect is to 
consider the turn on times as another renewal process driv- 
ing the spark discharge process. This type of model is 
known as a branching renewal process CRe-F. 30]. 

The noise processes that were used to develop this model 
were short term stationary in that the statistics o-F the 
process did not change over the observation interval . Not 
all sources o-F power noise are stationary even over a short 
observation interval. Physical effects such as wind, solar 
heating and varying line loads can act to make certain types 
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ot power line noise highly variable from observation to 
observation. Although these types of noise were observed 
during data collection, the data used to develop the model 
was taken -from noise sources that were stationary tor the 
length of the ten minute data records . 

C. OBSERVED TIME-DOMAIN CHARACTERISTICS 

The gap noise sources tor this study were observed on 
utility distribution lines in the vicinity ot the Naval 
Postgraduate School and were chosen to illustrate parameters 
ot the noise model. Fig. 8 is a typical time-domain observ- 
ation ot a gap discharge process observed at 3 MHz and 
envelope demodulated with a 10-kHz Gaussian bandpass tilter. 
The important character i st i cs to note are: 

1. the process has an on-ott-on modulation at a 120- 
Hz rate related to 60-Hz wavetorm ot the power 
line, and 

2. the pulse groups that result trom the modulation 
have a probabilistic number ot impulses occurring 
in each group and a probab i 1 i st i c interarrival 
time between pulses in a group. 

In the tirst pulse group, 9 pulses occur with varying 

amplitudes and interpulse arrival times. In the second 

pulse group starting approx i mate ly 8.33 msec later, 8 pulses 

occur, again with varying amplitudes and interarrival times. 

One group ot pulses is associated with either the positive 

or negative polarity ot the line voltage wavetorm and the 

other group with the opposite polarity. I dent i t i cat i on ot 
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Figure 8. 



Envelope oF Demodulated Gap Noise 
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the polarity is impossible without physically locating the 



source . 

The average amplitude o*P the impulses will be one o *P the 
parameters tor the model. An estimate of this parameter can 
be made from the above presentation. It is important to 
note that tor this highly impulsive type ot noise that the 
observed amplitude is a tunction ot the shape and width ot 
the ettective bandpass ti Iter and the detector char- 
acteristics. 

At this point it should be noted that the model can be 
approached at two levels ot complexity. In the simple 
version the data trom both the positive and negative phases 
is considered as a whole and averaged over the two phases to 
determine a set ot average parameters. In the more complex, 
but more accurate model , two* sets ot model parameters are 
determined - one tor each phase. For the remainder ot this 
work the more complex case will be considered. The results 
tor the simple case are presented in Moose and O'Dwyer CRet. 
31 1 and can easily be determined trom the complex case by 
setting the parameters equal tor both phases and simplitying 
the resulting equations. 

In order to better characterize the interpulse arrival 
times, which will be used to determine two model parameters, 
the rising raster capability described in the instru- 
mentation section was used to generate the display o*P a 
di-P-Perent gap noise source which is shown in Fig. 9. In 
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Figure 9 




Multiple Scans of Gap Noise Interference 
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•this picture the amplitude o-F each individual record was 
reduced and thresholded so that all that essentially remains 
is the time o-F impulse arrival data -For 28 time records. 
The time base o-F the display was intentionally synchronized 
to the power line to -Facilitate taking data and this 
accounts -For the regularity o-F the pulse groups -From 
observation to observation. In this view the interarrival 
times -For 56 pulse groups (28 groups o-F one polarity and 28 
o-F the opposite polarity) can be determined along with the 
number o-F pulses in each o-F the 56 successive pulse groups. 
The average number o-F pulses per group considered separately 
-For each phase will also be used as estimates o-F two model 
parameters . 

Fig. 10 is a histogram o-F the distribution o-F the 
interarrival times between the observed impulses -For the 
same noise process shown in Fig. 9. To use this data to 
generate parameters -For the noise model, the histogram will 
be approximated by a continuous density -Function. The gamma 
density -Function CRe-F. 32:p. 18] 

r - 1 - ( Xt ) 

p(t) = X(Xt) e / ( r- 1 ) ! (3-1) 
was chosen because it showed a close -Fit to the data, and 
its character i st i c -Function, which will be used in later 
derivations, was particularly simple. The character i st i c 
-Function o-F the gamma density is 

-r 

0(jo>) = (1 - j(uX) . (3-2) 
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Figure 10. Histogram of Impulse Interarrival Times (Case 1) 



The r and X parameters are simply interpreted in terms 
o-P an underlying Poisson process as the time to the rth 
point o-P a Poisson process o-P intensity X. I-P p is the 
random variable assigned to the observed interarrival times 
then the unbiased estimate o-P the mean o-P p is 

N 

est(p) = (1/N) E pi (3-3) 

i = l 

and the unbiased estimate o-P the variance is 

2 N 2 _ 2 

est ( <J M ) = -C E pi - N est(p) >/(N-l). (3-4) 

i = l 

The unknown parameters o-P the desired gamma density function 

are estimated using the method o-P moments CRe-P. 33 : p . 2503 

using the following relations; 

_ 2 2 2 
r = est(p) /est(a M ) and X = est(p) /est (a M ). (3-5) 

Using this approach, the gamma density function used to 

approximate the interarrival time histogram is also plotted 

on Fig. 10. Figs. 11 and 12 are based on the same data as 

Fig. 10. However, they are separated according to di-P-Perent 

phases. It can be seen -Prom the r and X parameters in Table 

III that there is a signi-Picant difference between phases o-P 

the same gap noise process shown in the histogram even 

though the mean value o-P each data set was almost the same. 

The data in Fig. 11 was more clustered than the Fig. 12 

data. This results in a higher r and lower X parameter for 
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Histogram of Interarrival Times (Case 1 - Phase A) 



GAMMA DENSITY FUNCTION, N=96 
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Figure 12. Histogram of Interarrival Times (Case 1 - Phase B) 



"the same mean value. 



Figs. 13 and 14 are additional sets oF 



gap noise data that were observed during the research. 

Table III summarizes the r and X parameters -Pitted to 
the data described above. 

r X 



Case 


i 






33.2 


37,000 


Case 


i 


Phase 


A 


63.4 


71 ,400 


Case 


i 


Phase 


B 


23.3 


25,900 


Case 


2 






44.3 


74,600 


Case 


3 






14.1 


19,200 



TABLE III 

A -Pinal model parameter will be called T & and is a constant 
delay or o-P-Pset o-P every second pulse group that exists with 
respect to the Fundamental voltage waveform. 

In summary the inputs to the model based on time-domain 
data are; 

1. An estimate of the average amplitude of all 
observed pulses. 

2. An estimate of the average number of pulses per 
group For each phase. 

3. An estimate oF the mean oF the interarrival time 
oF the pulses For each phase. 

4. An estimate oF the variance oF the interarrival 
time between pulses For each phase. 

5. An estimate oF the delay in the start time oF the 
pulse groups oF one phase with respect to the 
Fundamental voltage waveForm. 
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Figure 13. Histogram of Interarr i val Times (Case 2) 
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Figure 14. Histogram of Interarrival Times (Case 3) 



The average ampl ltude oF pulses -Prom diFFerent gap noise 
sources is highly variable, ranging -Prom “the instrumentation 
noise floor “to higher “than any observable signal in -the HF 
band. The average number o-P pulses per group has been 
observed -to vary -Prom 1 to greater than 20. The mean o-P the 
interarrival times ranges -Prom .1 to 1 ms and the 
coefficient o-P variation (standard deviation divided by the 
mean) o-P the interarrival times has ranged -Prom approx- 
imately .1 to .3 for the gap noise processes we have 
observed . 



D. MODEL DEVELOPMENT 

The principal goal of this chapter is to describe a 
noise model that, with suitable choice o-P parameters, pro- 
vides an adequate simulation o-P the actual physical noise 
mechanism o-P gap noise. In this section the time-domain 
data and parameters derived in the previous section will be 
integrated into the Filtered impulse model developed in 
Chapter II. Recalling From Chapter II, a Fi ltered impulse 
noise process is given in complex envelope notation by 

N ( t ) -jQt 

n(t) = g(t) + £ a t e h(t-ti ) (3-6) 

i = l 



where p(a) is the amplitude distribution oF the 
N(t) is a unit counting process whose statistics 
the impulse arrival times, 0 is a random variable 



i mpul ses , 
determ i ne 
uni Formly 
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distributed between 0 and 2n 



and hCt) is the impulse 



response of the receiver f 1 Iter . 

In order to physically justify a -Filtered impulse model, 
the impulse duration must be small compared to the inverse 
bandwidth o-P the receiver -Filter. This condition is easily 
met in the case o-P gap discharges. Laboratory analysis o-P 
temporal char acter i st 1 cs o-P gap discharges -For various 
geometries shows that the impulse durations range -Prom 10 to 
100 's o-P nano-seconds CRe-F. 34]. Therefore, -For -Filter 
bandwidths up to 1 MHz the output noise process will only be 
a -Function o-P the incident time o-F the impulse and the 
-Filter response , not the waveform o-F the impulse. 

The specification o-F the statistics o-F NCt), the count- 
ing process that drives the model is then crucial to 



obta ini ng 


an accurate representat l on o-F the physical noise 


process . 


Based on the electrical char acter l st l cs o-F a gap 


di scharge 


and the observed time-domain characteristics, the 


-Po 1 1 ow l ng 


assumptions are made; 


1 . 


A primary series o-F event times separated by an 
interval T 0 /2 exists. T 0 is the -Fundamental 

period o-F the power line voltage waveform. The 

distribution o-F the time to the -First primary 
event is uniform over 0 to T 0 - 


2. 


A subsidiary process commences at the primary 
event times i.e. T 0 /2 , T 0 , 3T 0 /2 ... The 

subsidiary process is a renewal process that 
continues -For N1 or N2 renewals alternating 
between N1 and N2 renewals at successive primary 
po l nts . 
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3. The interarrival times to the -First and subsequent 
points o-F the subsidiary renewal process -For each 
phase are all independent and identically 
distributed (IID) within a pulse. Each phase has 
its own set o-F parameters. 

4. The weighting or amplitude distribution o-F all 
impulses in the subsidiary process is a constant. 

5. A constant o-P-Pset or delay designated by T D may 
exist at every second primary renewal point. 

The -First assumption is supported by the periodicity in 
the data that is related to the -Fundamental -Frequency o-F the 
power line. The large majority o-F observed gap noise sources 
had pulse groups on both the positive and negative phases o-F 
the -Fundamental waveform. This accounts -For the T 0 /2 per- 

iodicity. 

The second assumption addresses the differences between 
the sparking phenomenon on the positive and negative phases 



o-F "the 


-Fundamental . 


It is -Frequently 


observed that 


the 


ECNp o . D 


i s 


d i fferent 


than ECN NCa ] where 


Np o s and N* e a 


are 


the random 


number o-F 


discharges in the 


pulse groups. 


Thi s 



can be due to two effects; an asymmetrical gap geometry and 
the -Fundamental physical difference between the sparking 
mechanisms -For positive to ground and negative to ground 
sparks. Two deterministic integer constants; N1 and N2 
nearest to ECNpos ] and ECN MK «J respectively are assumed as 
the approximation to the variable number o-F sparks per ha 1 -F 
cycle o-F the -Fundamental waveform. 
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The -third assumption is based on the -Fact that the 
"inception o-F gap discharges in natural air and the 
development o-F electron avalanche are -Fundamentally 
probabilistic processes that depend on atmospheric pressure, 
humidity, presence o-F natural ions, electrode surface and so 
on" CReF. 34D. In view o-F the above statement and 
considering the empirical data, the j ust i -Fi cat i on -For devel- 
oping the gap discharge as a probabilistic process is well 
-Founded. The assumption o-F independence -From discharge to 
discharge is not as well justiFied. EFFects such as 
electrode heating a-Fter the initial discharge in a cycle 
could act to make the average interarrival time vary -From 
pulse to pulse within a pulse group. However, to develop a 
tractable model the assumption o-F independence between 
arrival times is made. The gamma density -Function deter- 
mined by the estimate o-F the mean and variance oF the 
interarrival times For each phase is used to deFine the PDF 
oF the interarrival times. 

The Fourth assumption concerns the amplitude oF the 
impulses. In most cases the amplitude oF the impulses 
within a pulse was nearly constant. In some cases there was 
a diFFerence in the average amplitude From negative phase to 
positive phase, however to keep the model tractable the 
amplitude oF all impulses is assumed constant. 

The FiFth assumption incorporates an eFFect that will 
explain one oF the Features observed in the power spectral 



61 



density (PSD) o-P the observed process. The physical 

rationale behind this assumption is that there may exist a 
difference in the threshold between the positive and 

negative phases, thus consistently delaying or offsetting 
the start o-P the renewal process on one of the phases. 

The above assumptions place this model in a class o-P 
processes known as branching processes with the main process 
being a degenerate renewal process with an interarrival PDF 
o-P 

p(t) = <S (t-T 0 /2) (3-7) 

The subsidiary process is a renewal process of N1 or N2 

points. Consequently, the impulse source for the gap noise 

model is completely specified by the -following parameters: 
a - Amplitude o-P the impulses, 

N1,N2 - Number o-P impulses in the subsidiary process, 
p t - Mean value o-P impulse interarrival time, 

Var (p t ) , Var (p c ) - Variance o-P interarrival times, 

T D - O-P-Pset on one phase. 

To complete an interference scenario N 0 B C ff/2 and h(t) must 
also be speci-Pied. The expression for the branching renewal 
process is then 

M ( T ) N. -j8.„ 

n(t) = g(t) + E E a. „ e h(t-t.„-T.) (3-8) 

m= 1 n= 1 

Fig. 15a, b shows sample realizations o-P the envelope, 
E(t), o-P the noise process de-Pined above. 
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ECt) 




ECt) 




Figure 15. Envelope and Counting Function oF n<t) 



The Gaussian term was set to 0 and h(t) was approximated as 
a decaying exponential. Note that in comparison to an 
actual noise process illustrated in Fig. 8 the amplitude of 
the impulses is constant with a Fixed number of impulses on 
the positive and negative polarity of the line voltage 
waveform. Fig. 15c is the primary counting process, MCT), 
driving the subsidiary process shown in Fig. 15b. Fig. 15b 
includes a non-zero T D and its relationship to MCT) is shown 
in Fig . 15c . 

Fig. 16 shows a realization of the envelope of this 
process For N1 and N2 = 3, T = 30 ms and gamma density 
Function parameters rl and r2 = 32 and XI and \2 = 37,000. 
Random interarrival times with a gamma distribution were 
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Figure 16- Synthetic Gap Noise 



approximated by summing exponentially distributed random 



variables CReF. 35] . The impulse amplitudes were normalized 
to one and a background level oF quadrature Gaussian noise 
at NoB/2 = .0001 was added which simulates either receiver 
noise or high density, low amplitude impulsive noise. 

E. COMPARISON OF CALCULATED AND OBSERVED DATA 

One partial description of a noise process is the 
amplitude probability distribution (APD) of the envelope. 
This description has been used extensively in the analysis 
of noise and in deriving optimum receivers. For the NCt) 
specified in the previous section this calculation would be 
diFFicult. A second noise process descriptor is the 
spectrum of the noise process. This descriptor is suited to 
our model where the structure of the noise is contained in 
the counting process that drives the impulse generation. In 
addition to receiver noise performance evaluation, spectral 
analysis can also be used For noise source identification 
and isolation CReF. 36]. 

In order to determine the spectrum oF the envelope 
squared oF the branching renewal process deFined by Eqn. 
3-8, the spectrum oF the envelope squared oF a non- branch l ng 
renewal process is determined in Appendix D. The results 
obtained in Appendix D are then used to derive the spectrum 
oF the the branching renewal process which is postulated to 
model gap noise. Starting with the branching renewal 
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process given in Eqn. 3-8, and assuming g(t) is equal "to 
zero, it is shown in Appendix E -that -the average o-F the 
Bartlett estimate of the PSD is: 

2 5 

S ( a) , T ) = | Ha C (i) ) | /T L Term (n) (3-9) 

n= 1 

where n is an index to the 5 Terms defined below. These 
-Five Terms involved in the summation in Eqn. 3-9 are de-Fined 
as -Fol 1 ows : 

Term 1 - The inter-pulse group summations between pulses in 

pulse groups with N1 pulses 

Term 1 = M2 5 <(d,N 1) (3-9a) 

Term 2 - The inter-pulse group summations between pulses in 

pulse groups with N2 pulses 

Term 2 = M2 5 (cd,N 2) (3-9b) 

Term 3 - The intra-pulse group summations between pulses in 

groups with N1 pulses 

cos(M2 2nu)/(D 0 ) - 1 

Term 3 = < - M2 > n<ui,Nl) (3-9c) 

cos (2no)/(i) 0 ) - 1 

Term 4 - The intra-pulse group summations between pulses in 
groups with N2 pulses 

cos(M2 2n(i)/(j 0 ) - 1 

Term 4 = < - M2 > Q(cd,N 2) (3-9d) 

cos (2n(i)/(u 0 ) - 1 

Term 5 - The intra-pulse group summations between pulses in 

groups with N1 pulses per group to pulses in groups with N2 
pulses and vice versa 
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cos(M2 2n(i)/(i) 0 ) - 1 



Term 5 = 2cos ( ( tu ( T 0 /2-T D ) ) -C ■ 



■>^(0! , N 1 ,0 a ,N2) 



cos (2nu)/(i) 0 ) 



( 3-9e ) 

where 0 t and 0* are -the characteristic -Functions associated 

with the interarrival times -For each phase. The parameters 

N1 and N2 are the number o-F renewals associated with each 

phase and § , ft, and are de-Fined below: 

2 N 

N - N 1 0 ( j (i> ) | 0(jU)) (1 - 0(j(d) ) 

5(0, N) = - 2 Re-C >, 



|1 



2 

- 0( JO)) | 



2 

(1 - 0(jCD)) (3-10) 



and 



« (0 ,N) 



N+l 2 

|0(jO)) - 0(j<u) | 



2 

I 1 - 0( jw) | 



(3-11) 



U/(0, ,N1 ,0 a ,N2) 



Nl+1 N2+1 

(0i (j<d) - 0i (jtu) ) (0a (“JO)) - 0 a (-jd)) ) 



( 1-01 ( JOJ) ) ( 1 - 0a ("JW) ) 

(3-12) 



As a test o-F the model the PSD o*F the gap noise process 
shown in Fig. 9 was determined using a spectrum analyzer and 
also analytically determined using only the model parameters 



From the 


time-domain data 


and the methodology 


o-F 


the 


previous 


section. 


Fig. 17 


is a computer plot 


o-F 


the 


analytic 


est i mate 


of the PSD 


which is compared i 


with 


the 



observed PSD computed on a Wavetek UA500A spectrum analyzer 
shown in Fig. 18. 
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AMPLITUDE 



Figure 18 . 




Observed Estimate o-P Power Spectral Density (PSD) 
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Figure 17. Predicted Power Spectral Density (PSD) 

Comparing the two PSD's, it is seen that the analytical 
expression correctly predicted all of the significant 
features seen in the observed PSD. 

The components of the spectrum due to the different 
terms is of interest and offers some insight into the 
origins of the features of the observed spectrum. The dom- 
inant feature of the observed spectrum are the periodic 120 
Hz spectral lines. At 800 Hz the spectral lines change, and 
appear as odd harmonics of 60 Hz, however, still remaining 
at 120 Hz intervals. Terms 3-5 are responsible for the 
these spectral line features. Each of these terms is a comb 
function, determined by the fundamental frequency and 
observation interval, multiplied by an envelope determined 
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by the interarrival time distribution. Figs. 19 and 20 are 
Terms 3 and 4 - respectively and clearly show the phase to 
phase difference. The more peaked behavior oF Fig. 19 at 
1100 Hz is due directly to the more clustered interarrival 
time behavior on one phase. This effect was observed in the 
time domain in Fig. 11. When the terms shown in Figs. 19 
and 20 are added together, they Form 60 Hz harmonic spectral 
lines and when Term 5, shown in Fig. 21, is added, positive 
reinForcement occurs at 120 Hz and negative at 60 Hz 
harmonics. This produces the observed 120 Hz harmonic 
spectral lines. When the cancellation is not complete due 
to diFFerences in the character i st i c Function t From phase to 
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Figure 21. PSD Term due to Phase A to Phase B Intra-Pulse Group 
Corre 1 at l ons 
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Figure 19. PSD Term due to Phase A Intra-Pulse Group 
Correlat ions 
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Figure 20. PSD Term due to Phase B Intra-Pulse Group 
Cor re 1 at i ons 
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phase, small 60 Hz components can be present. This effect 
is seen in both -the observed and predicted spectra. 

The observed spectra changeover -Prom 120 Hz harmonics to 
odd 60 Hz harmonics at 800 Hz is also correctly predicted by 
the model. This e-P-Pect is due to the factor, T D , appearing 



in Term 5 


which causes the 


odd 


60 


Hz 


harmonics 


to 


be 


positively 


reinForced and 


the 


120 


Hz 


harmon i cs 


to 


be 



canceled out above 800 Hz. The -Pinal effect shown in the 
model and seen in the observed spectra is the presence o-P 
continuous spectral components. The continuous spectra due 
to Terms 1 and 2 is shown in Figs. 22 and 23 and also 
differs From phase to phase*. The continuous spectra peaks 
markedly at 1100 Hz and this effect is apparent in the 
observed spectra as a rise in the noise Floor. The con- 
tinuous spectral term arises when Term 1 is added to Term 3 
and Term 2 to Term 4. The peak value at 1100 Hz in Terms 1 
and 2 is greater than the negative values in Terms 3 and 4 
at the same Frequency and when added together causes a 
continuous rise in the spectral Floor. The reason For this 
is that the i to i -+• 1 and i-1 interarrival times within a 
pulse group are more correlated than the correspond i ng 
interarrival times between diFFerent pulse groups. The more 
correlated arrivals appear in Terms 1 and 2 and cause the 
spectral peak to be larger. 
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Figure 22. PSD Term due *to Phase A Inter-Pulse Group 
Correlations 



o 




Figure 23. PSD Term due -to Phase B Inter-Pulse Group 
Correlations 
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F. 



SUMMARY 



A nine parameter model For a single source of gap noise 
interference was developed and shown to accurately predict; 
the PSD of narrow bandwidth envelope demodulated gap noise 
interference. IF less accuracy is desired, phase to phase 
diFFerences can be neglected and the resulting model has 
only Five parameters. The model was speciFied in terms oF 
an impulse driving Function, and techniques to estimate the 
model parameters From time domain data were described. 
Since the model speciFies a complex driving impulse Function 
it can be used as a simulated interFerence source For 
arbitrary Filters and receivers. 
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IV. CORONA NOISE MODEL 



A. INTRODUCTION 

In addition to gap noise which was discussed in Chapter 
III, another major source o-F power line radio noise is 
corona noise. Corona is di-F-Ferent in many respects -From gap 
noise, and the model used to describe it will take a 
di-F-Ferent approach. In contrast to gap noise, which is 
observed on both electric power transmission and dis- 
tribution lines, corona noise is observed only on power 
transmission lines and generally has a -Fundamental -Frequency 
o-F 180 Hz correspond i ng to the three phases o-F a 60 Hz 
alternating curr.ent (AC) system. Corona noise is caused by 
a partial breakdown o-F air surrounding a conductor which is 
at a high potential. Consequently, the impulsive structure 
o-F corona noise is not as well de-Fined as the highly 
impulsive structure o-F gap noise. Corona noise -From AC 
sources appears somewhat like amplitude modulated Gaussian 
noi se . 



B. CORONA NOISE THEORY 

When an increasing potential difference is applied 
between two electrodes, a breakdown voltage is reached that 
is character i zed by the transition o-F air -From a poor 
electrical conductor to a good conductor. I-F the -Field in 
the gap between the two electrodes exceeds the electric 
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strength of air, which is about 30 kV/cm -For gaps over 1 cm 
in length, prior to the spark breakdown voltage, then corona 
will exist. Sharp points and small radius wires are avoided 
whenever possible on power lines, and corona generally does 
not become a significant source of radio interference until 
voltages exceed about 70 kV CRef. 28 : p . 84]. 

As the voltage on an AC line is increased, corona will 
appear and be associated with either the positive, negative 
or possibly both half-cycles of the fundamental waveform. 
For an aged, clean transmission line it appears that the 
threshold voltage for negative corona is less than that of 
positive corona. If the transmission line has any unusual 
characterist ics , the threshold voltages of positive and 
negative corona become difficult to categorize CRef. 37]. 
In particular, it has been observed that positive streamer 
corona , which is more disruptive of communications fre- 
quencies than negative corona, predominates on newly 
installed lines and during periods of prec i p i tat i on . For 
modeling purposes, the important character i st i cs of corona 
noise are not a precise description of the physical 
mechanisms but the relationship between the observed inter- 
ference waveform and the line voltage. 

The procedure most commonly used for the prediction of 
interference due to corona noise gives the noise power as a 
function of the char acter i st i cs of the interference 
scenario. The random properties of the corona noise are 
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averaged out CRef. 383 . These characteristics would 
commonly include distance -From the line, “the line voltage 
and a weather correction term. Each of these terms is 
exper imental ly verified and related to some physical process 
in the generation of corona noise. While this procedure 
gives an estimate of the noise power in a specified 
bandwidth, it does not incorporate the probabilistic nature 
or time-domain periodicity observed in the noise. 

The approach taken in this chapter will continue with 
the filtered impulse theme of the dissertation. For co)ronai 
noise, the counting process , N(T), will be modeled as a 
periodically modulated stream of impulses. This approach 
will have the advantage of being accurately able to simulate 
the time- and frequency-domain behavior of corona noise. 
However, as with the gap noise model it will suffer from the 
disadvantage that for each scenario, a set of parameters 
must be experimentally determined. 

C. TIME DOMAIN CHARACTERISTICS 

The time-domain data for corona noise was obtained using 
the instrumented van described in Appendix A. The antenna 
was a 3-meter whip mounted on the van (which was parked 
underneath a 3-phase 500-kV transmission line). The weather 
conditions at the time of measurement were clear and sunny. 
Records of the corona noise were obtained using the envelope 
detected output of the HP-141T described in Appendix A and 
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recording the data on a 20-kHz bandwidth, 85-dB dynamic 
range digital audio system. The recording process allowed 
subsequent analysis in the laboratory. 

A 37-ms portion o-F a data record is shown in Fig. 24. 
The data was taken at a 300-kHz center -Frequency and 10-kHz 
intermediate -Frequency (IF). This particular picture was 
obtained by playing back the tape, re-digitizing the data 
and displaying it with a signal-analysis software package. 
One of the important -Features o-F the time-domain record in 
Fig. 24 is the periodicity. The -Fundamental -Frequency of 
180 Hz can be determined -From a periodic rise in the mean 
value o-F the noise and -From a periodic increase in high- 
amplitude impulses. The -Fundamental noise -Frequency o-F 180 
Hz on a three phase line indicates that the dominate noise 
is generated on one half-cycle o-F the -Fundamental 60-Hz 
waveform on each o-F the three phases o-F the line. This is 
consistent with low-level negative corona predicted -For aged 
lines in clear weather. 

The periodogram spectrum o-F the data record described 
above is shown in Fig. 25. The time-domain observations o-F 
the 180-Hz periodicity are confirmed and a 360-Hz harmonic 
o-F the -Fundamental is also present. At higher frequencies, 
no harmonics are observed and the spectrum becomes white. 
This is in distinct contrast to the gap noise described in 
Chapter III where the spectrum o-F gap noise has harmonics 
well into the kHz range. The periodicity can also be seen 
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Figure 24. Corona Noise xn Time Domain 
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Figure 25. Corona Noise in Frequency Domain 



in Fig. 26 (which is “the same data recorded on the 3-Axis 
display). The contribution oF each phase oF the power line 
to the total corona noise process can be clearly seen in the 
lower picture in Fig. 26. 

Another important Feature of the data For modeling 
purposes is that the receiver response to an individual 
impulse overlaps with the response oF preceding impulses. 
Although the noise is still impulsive, the simpliFying 
assumption oF non-overlapping pulses is not valid. It is 
seen in Fig. 24 that the amplitude oF the impulses is random 
and there is no Fine-grain interarrival time structure other 
than the periodicity discussed above. Both oF these 
observations are in direct contrast to gap noise where the 
amplitude From impulse to impulse is almost constant, and 
there is a deFinite impulse interarrival time structure. 
These observations provide a method to distinguish between 
gap noise and corona noise. 

An interesting aspect oF noise From transmission lines 
was noted while taking measurements For the corona model. 
It was diFFicult to obtain a recording oF power-line noise 
in the vicinity oF high voltage transmission lines due only 
to corona noise. In many cases, a gap noise component was 
superimposed on the corona noise. This can be seen by 
comparing Figs. 26 and 27. Fig. 27 is a view oF both gap 
and corona noise. The corona noise amplitude peaks at the 2 
pV level, and the amplitude oF the gap-noise impulses is 
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F 1 gure 26 . 



Corona Noise 



as shown on 3-Axis Display 
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Figure 27. Gap and Corona Noise 
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about 3.5 pV. The two lower pictures are ot identical data 



with the threshold varied to emphasize each noise type. The 
lower picture has the threshold set at about 3 pV 
eliminating the 180-Hz corona noise and showing the gap 
noise process occurring at a 120-Hz rate. This corresponds 
to a gap noise source sparking on the positive and negative 
halt-cycle ot a single phase ot a 3 phase line. When the 
threshold is lowered to about 1 pV, the 180-Hz corona noise 
dominates as seen in the middle picture ot Fig. 27. This 
phenomenon is due to corona discharge on the negative halt- 
cycle ot each phase ot a 3-phase system. Without the abil- 
ity to examine the time- and trequency-doma i n behavior ot 
the noise using a scanning analyzer and 3-Axis display, it 
would be very ditticult to d l tterent i ate between gap and 
corona noise and make accurate measurements ot either. 

Fig. 28 illustrates corona noise with ditterent phase' 
to-phase character i st i cs . In this example the noise power 
due to the corona is unequal trom phase to phase. The noise 
occurring on the dominant phase is approx imate ly twice the 
peak voltage ot the phase with the lowest amplitude. 

D. CORONA NOISE MODEL 

Based on the above observations ot the demodulated 
envelope ot corona noise, a specitic tiltered impulse model 
based on Eqn. 2-18 will be developed. This model is similar 
to models used in synthetic hydrology to account tor 
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Figure 28. 



Corona Noise 



Uneven Phase 



to Phase 
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seasonal variations in streamf low data CReF. 391- Assuming 
no Gaussian noise component, a lowpass complex -Filtered 
noise process can be written 



N(t) -je t 

n ( t ) = E a 4 e 6 (t-t* ) . (4-1) 

i = 1 



Since corona noise is a high, density process, the 
cannot be simply deFined as in Eqn. 2-22 since the 
oF the Fi Iter to adjacent pulses will overlap, 
cumvent this diFFiculty it will be assumed that the 
oF n(t) can be modeled by 



envelope 
response 
To cir- 
enve lope 



N(t) 

e ( t ) = E a t g(t-t 4 ) (4-2) 

i = l 



where both e(t) and g(t) are real, positive and g(t) depends 

on h(t) and the process intensity. 

The assumptions to de-Fine the model are as -Follows; 

1. The rate parameter -For the counting -Function N(t) 
is periodically modulated at -Frequency (do and will 
be approximated by a -Finite number o-F terms in a 
Fourier series. 



2. The amplitude distribution of successive pulses 
given by the a t 's are independent and identically 
distributed with density Function p(a). 

3. The start time of the periodic modulation is 

uniformly distributed over the period T 0 where T 0 
= 1 / (do - 



Assumption one is motivated by the periodicity observed in 
the data. In most cases the Fundamental period, u) 0 will be 
180 Hz. IF the noise is more dominant on one phase oF the 
line, as shown in Fig. 28, then 60 Hz components oF the 
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noise will exist and for a more accurate model the 
fundamental period would be 60 Hz. The rate parameter is 
given by 



N t 

\(t) = E a„ cos (u)o nt ) 
n=0 



(4-3) 



where N T is the order at which the Fourier series is 
truncated . 

The second assumption is arbitrary in terms of the 
distribution function chosen for the impulses. It is 
obvious that random impulse amplitudes are being generated 
by the cumulative effect of the corona noise sources. 
However, it would be difficult to accurately fit a density 
function to the data due to the time varying nature of the 
process. The independence of the amplitude bursts is also 
conjecture but facilitates the development of a tractable 
model. The Rayleigh density function is 



2 

-a /2fl 

p ( a ) = (a/6 ) e a > 0 

=0 a < 0 (4-4) 

and was picked to model the random amplitudes. It is 

analytically tractable and can be justified by assuming the 
quadrature components of the impulse amplitudes are in- 
dependent Gaussian events. 

The final assumption is made because in most 
interference scenarios, the power line noise process will 
not be synchronized to the process being interfered with. 
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Thus the phase will be effect 1 ve 1 y randomized over one 



-fundamental period. In some o-F the observations o-P corona 
noise Csee Figs. 26-28) this assumption was clearly violated 
and the 1 nst r umentat 1 on was 1 nt ent i ona 1 1 y synchronized to 
the power line For clarity oF presentat i on . 

The model could be Further generalized by allowing the 
amplitude density Function p(a) to be a periodic Function oF 
time, p(a,t) independent oF the rate parameter. This added 
complexity would make the model more accurate; however, more 
parameters would have to be estimated. 

The above assumptions classiFy this process as a 
compound non- homogeneous Filtered Poisson process. A con- 
venient tool to study these processes is the cumulant 
generating Function CReF. 41:p. 117] 

se (t ) 

14 / ( s) = ln-CECe ]> (4-5) 
where In is the natural logarithm. The joint cumulant gen- 
erating Function is 

S, e(ti ) + Sa e ( ta ) 

\|/(si , s 2 ) = ln< ECe ]>. (4-6) 



For the noise process deFined by the First two assumptions, 
it is shown in Appendix F that 



\\j ( s ) 



00 



p ( a ) 



00 

X ( a ) 



sah ( t-a ) 

C e 



1 ] dada 



and 



(4-7) 
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\j/ ( S| , Sa ) 



a> 



oo 



p(a) 



X C a ) 



u u 

"00 “00 



Si ah (*ti -a ) 

C e 



s a ah ( *be -a ) 

- 1 3dada - 

(4-8) 



From the properties oF the cumulant generating Functions it 
can be shown that 

ECe(t) 3 = y' (0) , (4-9) 

Var C e ( t ) 3 = i|/''(0) (4-10) 

and 

2 

c ly ( 0 , 0 ) 

Cov(t t ,t a ) = . (4-11) 

DSi *OS a 

It is straightforward to show that 



ECe(t) 3 



E C a 3 



00 

1 N t 

E a„ cos (u) 0 na)h(t-a)da , 

J n= 1 



(4-12) 



and 



Var Ce (t ) 3 



2 

ECa 3 



00 

1 Nt 2 

E a n cos ( (Jo na ) h (t-a)da 

J n= 1 



-00 



(4-13) 



CovCti ,t a 3 



2 

ECa 3 



00 

1 Nt 

E a„ cos ( (Jo na ) h ( ti -a ) h ( t 2 - a ) da . 

J n= 1 



- ao 

(4-14) 



The autocorrelation Function can be Found in terms oF the 
above Functions as 



R(ti,t. ) = Cov(ti ,t a ) + ECe(ti ) 3ECe(t a ) 3 . (4-15) 
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The above equations, prior to -the application o-P 
assumption 3, de-Pine a wide sense cyclostat ionary process 
CRe-P. 40] where 

ECe(t + mT 0 ) ] = ECe(t)] C4-1S) 

and 

RCtj + mT 0 ,t a + mT 0 ) = R(t t ,t a ) . (4-17) 

The statistics o-P the process are invariant to a shi-Pt o-P 
the time origin by integral multiples o-P the period, T 0 . 
This character i zat ion is particularly appropriate when the 
process is observed in synchron i zat i on with the fundamental 
power line frequency. 

When assumption 3 is applied, the process start times 
are randomly shi-Pted over one period, the process becomes 
wide sense stationary and the mean and autocorre 1 at i on 
become 

To 

ECn (t ) ] dt 



ECn(t)] = 



(4-18) 



0 



and 



R(t) = 



To 

R ( t+T , t ) dt . 



(4-19) 



0 



The power spectral density (PSD) o-P the process is 



S ( a) ) = 



_ “ JUIT 

R (t ) e dT . 



(4-20) 
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The -Following algorithm is an example o-F how -the noise 
process speci-Fied above can be digitally generated. For the 
algorithm speci-Fied below; N 0 is the number o-F samples in 
one period To , <n t > is the real sequence to be generated, i 
is the index variable -For the sequence and j and k are dummy 
index variables. 

1 . i * 0 

2. n ( i ) = 0 

3. Generate a Rayleigh random variable, R, with 
parameter 8 and set n(i) = n(i) + R 

Ny 

4. Set intensity X = E a„ cos(o) 0 ni/N 0 ) 

n=0 

5. Generate a Poisson random variable, P, and set j = 
trunc(P/X + .5) 

6. Set n(k) = 0 -For k = i + 1 to j 

7 • i “ j 

8. Go to Step 3 

The above algorithm will continuously generate a modulated 
Poisson impulse process that may be digitally -Fi ltered 
concurrent with the impulse process generation or subsequent 
to the generation o-F a complete time record . 

Two additional points should be noted. A number o-F 
impulse amplitudes could stack up on each other i-F the 
truncation operation returns a zero value. This corresponds 
to the physical case where multiple impulses occur and are 
unresolved within one sampling period. Before a record o-F 
noise samples generated by the above process is used a 
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truncated uniform random variable between 0 and N 0 , the sampled 

period, should be subtracted from i. This ensures “that: the 
start time is uniformly distributed over one period. 

E. PARAMETER ESTIMATION 

Having specified the model, the task of determining the 
parameters for a given interference scenario remains. The 
parameters to be determined are the *Ca„>, the Fourier Series 
coefficients of the periodically varying rate function and 
the A parameter in the probability density function of 
impulse amplitudes. As an example, the parameters for the 
noise process described in the time domain observations 
Section will be determined. 

One method of estimating these parameters is by 
considering the predicted mean and variance of the process 
as given by Eqns. 4-12 and 4-13 and matching these 
parameters to the data. In order to do this analysis, the 
corona source has to be observed as a cyclostationary 
process. This requires synchronizing the data collection 
instrumentation to the line frequency so that the same point 
in each cycle can be examined. This is difficult to do for 
field measurements . 

Another option for determining the model parameters is 
the predicted stationary mean, variance and PSD or auto- 
correlation function. The PSD as defined in Eqn. 4-20 is 
stationary so no line synchronization is required and 
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instrumentation is readily available to make a periodogram 

estimate oF the PSD. 

For the measurements taken in this research, the lowpass 
impulse response of the Filter was approximately 

a 

-at 

h(t ) = (a/n )e t 2 0 

= 0 t < 0. (4-21) 

With the above definition the Filter has unity gain at u) 

equal to 0. The eFFective Fi Iter response g(t) was approx- 

i mated by h (t ) . 

In Appendix G it is shown that 

2 2 
-u)o n 



1/2 N t 4a 

ECe(t)!3 * (Bn/2) E e a„cos(u) 0 nt) (4-22) 

n=0 

1/2 

where (nfl/ 2) is the First moment oF p(a), the probability 

density Function oF the impulse amplitudes. The stationary 
mean is 



1/2 

E C e ( t ) 1 * (Bn/2) a 6 . (4-23) 

The covariance oF e(t) is shown to be 



Cov (tt , t* ) 



2 

-a ( ti -t # ) /2 

(a/n ) 2Be 



CD 2 

' N T -2au 

E a„cos(a) 0 n(u+(ti +t, )/2)e du . 

J n=0 

1 1, -t, | (4-24) 
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Using Eqns. 4-15 and 4-19 and substituting the above 
results, the stationary autocorre 1 at i on -Function is 



2 

-oct / 2 1/2 2 

R(t) = 20 Ce a 0 (oc/8n) + (l/4)a 0 



2 2 
- (J 0 n 



N t 

+ (1/8) E e 
n= 1 



4a 2 

a 0 cos ( (i) 0 nT ) ] . 

(4-25) 



Taking the Fourier transform, the resulting PSD is 



S ( u) ) 



2 

1/2 -w /2a 

20 C(2) a 0 e 



2 

(n / 2 ) a 0 6(u) 



2 2 
- (D 0 n 



Nt 

+ ( n /8 ) E e 
n= 1 



4a 2 

a„ <5 ( w ± n(u 0 ) . 

(4-26) 



The terms in the PSD in order are; (1) a white noise term 
multiplied by the Frequency response oF the Gaussian Filter, 
(2) an impulse at tu = 0 From the mean value oF the process 
and (3) periodic components at harmonics oF 180 Hz due to 
the modulation oF the driving Function. 

To estimate the model parameters, the sample mean and 
mean square estimates were used in conjunction with Eqns. 4- 
23 and 4-24 to solve For a 0 and 0. The penodogram estimate 
oF the spectrum was used to determine the magnitude oF the 
a t and a a co-eFFi c i ents . The Fourier series was truncated 
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at N t 



= 2 since no higher order -terms were observed in Fig. 



25. For -the data taken in Fig. 24 and 25 -the parameter a in 
Eqn . 4-21 was 7.11xl0 9 . Table IV lists the parameters of the 
assumed model. 



Coefficient of Gaussian Filter a 7 . 1 1 x 1 0 9 

Amplitude Density Parameter G 1.58xl0 21 

Fourier Co-efficients of Rate a 0 22,300 

Funct i on 

a t 5 , 700 

a a 2,500 



TABLE IV 

These parameters were inserted into the algorithm described 
previously For generating corona noise, and Five cycles oF 
the synthetic noise are shown in Fig. 29. Comparing the 
actua 1 noise in Fig. 24 and the generated noise in Fig. 29, 
the generated noise correctly models the random, impulsive 
character oF the actual noise and also its periodicity. The 
synthetic noise, however, does not incor-porate the constant 
bias above the 0 level seen in the actual noise. 

F. SUMMARY 

This chapter has presented a detailed analysis oF the 
time- and Frequency -doma i n behavior oF corona noise. A 
model For synthetic corona noise was postulated based on a 
Filtered non- homogeneous Poisson process. A methodology For 
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Figure 29. Synthetic Corona 



estimating the model parameters was developed. The ad- 
vantage oF this model over previous work is that it allows 
an accurate representat i on of the time-domain behavior of 
corona. The disadvantage of this model is that it is 
scenario dependent and For each interference situation a new 
set of model parameters must be specified. One limitation 
was introduced into the model by the assumption leading to 
Eqn. 4-2. The lowpass equivalent Filter is a Function oF 
the bandpass Filter and the process intensity, and 
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consequent ly , 
■the gap noise 



the model is not as general or as accurate as 
mode 1 . 
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V. ENERGY DETECTION R ECEIVERS 

A. INTRODUCTION 

Energy detection receivers in the high frequency CHF) 
radio band operate in an impulsive radio noise environment. 
To predict and improve the performance of energy detection 
receivers, accurate models o-P the signal and noise envir- 
onment must be used. The noise models developed in earlier 
Chapters will be used to simulate atmospheric and man-made 
radio noise en v 1 ronment s . Two types o-P energy detection re- 
ceivers will be considered: a -Past Fourier transform (FFT) 
processing system and a compressive receiver. To improve 
the performance o-P these receivers the concept o-P the 
locally optimum receiver is introduced to suggest a robust 
post -detect ion processor . 

B. ENERGY DETECTION 

The energy detector receiver measures the energy in a 
signal over a specific time interval. For this chapter we 
will consider the signal to be a modulated sinusoid, o-P un- 
known Frequency at the receiver. The signal is then 
represented as 

s(t) = Acos(iD a t + 9(t)) (5-1) 
where A and (u 0 are random variables and 0(t) is a slowly 
varying function so that the signal power is effectively 
confined to a narrow bandwidth B. 
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Defining the problem as a hypothesis test 
H 0 : xCt) = n(t) noise alone 

H ! : x(t) = n(t) + s(t) signal plus noise. (5-2) 
The well known energy detection receiver [Ref. 42] uses t e D , 
the received energy as the test statistic such that 



Under H c : t E o 



T 

1 2 

n ( t ) dt 



0 



Under H, : t e 0 



T 

■v 

Cn(t) 

V, 

0 



2 

+ s (t ) ] dt . 



(5-2) 



For a sampled bandpass process the test statistic is 



N 2 2 

H 0 : t e o = E C n c » + n 3 i ] 

i = l 



and 



N 2 2 

Hi : t e o = E C n c t + 5 ct ] + [n 3i + s s i 1 

i = l 



(5-3) 



x C i = n c i + s ci or x C i = n c i 
and likewise for the quadrature component. 

Fig. 30 shows the block diagram of a receiver that cal- 
culates this statistic. This receiver will be called a 
square and sum receiver. The square and sum receiver is the 
optimum receiver for detecting random phase signals in 
Gaussian noise at small signal levels. At large signal to 
noise ratios (SNR) this detector is very close to optimum 
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Figure 30. Square and Sum Receiver 

[Ref. 9:p. 370]. When N = BT and the signal is sampled at 
the Nyquist rate, the test statistic is the energy in the 
signal in a period T. 

When it is desired to implement the receiver shown in 
Fig. 30 over a wide frequency band relative to the target 
signal modulation bandwidth B, a number of options may be 
pursued. The channelized receiver approach in which an 
independent receiver operates at each frequency will often 
be too expensive. Two other possible methods for generating 
the test statistic are; (1) a compressive receiver as a 
spectrum analyzer CRef. 43] and (2) a fast Fourier transform 
CFFT) signal processor. 
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It has been shown that an FFT signal processor exactly 
implements the test statistic given in Eqn. 5-3 tor the 
assumed signal in a single FFT filter [Ref. 44]. Whether 
successive samples will be independent or not depends on the 
number of transforms per second relative to the signal 
bandwidth B. If the sample interval is greater than 1/B 
then successive samples will be independent. 

In the compressive receiver, shown in Fig. 31, a fast 
sweeping local oscillator linearly scans the band of 
interest. Any narrow bandwidth signal in the band will 
appear at the mixer output as a linear, frequency modulated 
CFM) , chirp signal. This chirp signal is passed through a 
weighting filter to reduce sidelobes and into a dispersive 
delay line (DDL) which is matched to the inverse of the 
linear FM sweep. The output of the DDL is envelope detected 
and, when referred to the sweep time, it provides an 
estimate of the spectrum of the input signal. It is easily 
shown that for a sinusoid, the detector output is 
proportional to the signal magnitude. Thus, the compressive 
receiver approximates the square root of the ith value of 
the test statistic in Eqn. 5-3. The samples will be 
independent if the scan revisit time is greater than 1/B. 

C. SQUARE AND SUN RECEIVER 

In this Section the energy detection performance of the 
square and sum receiver, shown in Fig. 30, will be examined 
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Figure 31- Block Diagram of Compressive Receiver 



as the noise statistics depart from a Gaussian distribution. 



The square and sum receiver, which is optimum for small 
signals in Gaussian noise, will be used as the reference 
receiver so that the performance of improved receiver 
structures may be evaluated. 

The number of samples, N, chosen for this study was 10. 
This is a small sample size and, for impulsive noise, the 
resulting test statistic will have a distinctly non-Gaussian 
distribution. In Fig. 32 the probability of false alarm 
(P FA ) is plotted for unit root -mean -square (RMS) Gaussian 
noise and the Hall models with V D equal to 2.1, 3, 5 and 12. 
(see Table II). The normalized threshold is in reference to 
the unit RMS power of the noise. 

A Monte Carlo simulation technique was used to obtain 
the data for the plot. Since the data are independent 
Bernoulli trials with a large number of trials, the norma 1 
approx i mat i on to the binomial was used to determine the 
confidence interval. For estimates of P P A greater than 5 x 
10‘ 3 with 20,000 trials, the true value of P P A will be 
within 20 percent of the estimate with a confidence of 95 
percent CRef. 55 : p . 282]. The noise samples were generated 
using the inverse method [Ref. 45 : p . 951] and the inverse 
cumulative distribution functions (CDFs) in Table I. 

In Fig. 33 the probability of detection ( P 0 ) curves are 
plotted for P P A equal to 1x10" a . A rejection method [Ref. 
45 : p . 952] was used to generate the random deviates since 
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Figure 32. P* a For Square and Sum Receiver in Gauss and Hall 
No i se 
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•the CDFs oF the signal plus noise densities were not Found. 



The P t> is plotted as a Function oF SNR For the same 
parametric conditions as used in Fig. 32. As anticipated 
the P D decreases dramatically as the noise becomes more 
impulsive. In the case where V D equals 12, the P D decreases 
From greater than 99 percent to 1 percent at a 5 dB SNR For 
the same P P * . 

A receiver that is optimum in Gaussian noise suFFers a 
major degradation when the noise statistics become non- 
Gaussian. One possible approach to improving the per- 
Formance is to design a receiver that is based on the 
statistical parameters oF the noise. A parametric approach 
will be discussed in the next section. 

D. LOCALLY OPTIMUM RECEIVER 

The concept oF a locally optimum receiver provides a 
general methodology For nonlinear receiver design in the 
presence oF non-Gaussian noise. The approach is to 
determine a receiver structure that is optimized For the 
diFFicult small signal case and then examine, and lF 
necessary, modi Fy the structure For the large signal case. 
In this Section the locally optimum receiver For energy 
detection in the presence oF Hall type impulsive noise will 
be determined under the assumption oF weak signal levels. 
Previous work has considered known signal detection in 
Cauchy noise which is related to one oF the Hall models 
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CRe-P. 463. 



After optimizing the design -for weak signal 



levels the performance o-P the receiver will be examined -for 



will be parametric; that is it will depend on the noise 
density -function parameters. 

To illustrate the method -for deriving the locally 
optimum receiver, consider the -following binary hypotheses 
test. We are given a sequence, x 4 , o-P N samples. Under the 
null hypotheses H 0 , we assume that the sequence consists 
only o-P independent and identically distributed noise 
observations with a common density p(n). Under the 
alternative hypothesis Hi , we assume that the sequence 
consists o-P an additive mixture o-P a small signal and noise. 
The problem may be stated as 
H 0 : Xi = n 4 noise alone 

H t : x 4 = n 4 + Ast signal plus noise 1 = 1 , . . ,N 

with A an unknown small amplitude. 

The optimum receiver, in terms o-P maximizing the 
probability o-P detection while keeping the probability o-P 
-false alarm below a certain level is given by the Neyman- 
Pearson lemma. This receiver compares the likelihood ratio 
test statistic to a threshold which is chosen to achieve the 
desired probability o-P -false alarm. 

The liklihood ratio test is 



moderate to large signal levels. The receiver structure 



Po < x | Ho ) 



> T dec 1 de H 0 , or 



A(x ; A) 



(5-4) 



Pi < x I Hi ) 



< T decide H t 
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where p 0 (x| H 0 ) and p t (x| Ht ) are the probabi 1 lty density 
■Functions -For the observation vector x under H 0 and H t 
respect 1 ve ly - Since the noise samples are independent and 
identical ly distr i buted , p a (x | H 0 ) and p t (x| H t ) can be 
written as products o-F the univariate noise densities: 

N 

p 0 (x) = n pt ( x i ) (5-5) 

i = l 



and 



N 

Pi (x) = n pi (x t - A s i ) . ( 5 - S ) 

i = l 

The test can then be written 



A(x; A) 



N 

n Pi (x* - 

i = l 



N 

n Pi (x t 

1 = 1 



Asi ) 



> T dec l de H a , or 



< T decide H t . (5-7) 



In deriving the locally optimum receiver several 
approaches have been used. Rudnick CRe-F. 473 derived a -Form 
o-F the locally optimum receiver by expanding the likelihood 
ratio in a Taylor series about x and truncating after the 
second order terms. Capon CRe-F. 483 showed that the locally 
optimum detector maximizes the slope o-F probability o-F 
detection versus signal strength -Function at the origin and 
employs a statistic based on the -Following test 

•o A ( x ; A ) 

t l o ( x ) = (5-8) 

dA A => 0 . 
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In Appendix H it is shown that the locally optimum 



receiver tor the detection of a random sinusoid is 

N 

t ' uo (r) = E g ( r i ) (5-9) 

i = l 

where g(*) is the zero-memory nonlinearity (ZMNL) defined 
be 1 ow 



f ' ' (r) f ' (r) 

g ( r ) = + (5-10) 

f ( r ) r f ( r ) 

and r is a vector of envelope samples. The function, f(r), 
for the Hall model of atmospheric noise discussed in Chapter 
I I is 



0-1 



f (r) = 



* Ky 

2 2 (0 - 1 ) /2 
( r + y ) 



(5-1 1 ) 



where K is a constant. From Eqn. 5-10 the ZMNL of the test 
statistic is 

2 

(0 + 1 ) ( 0 + 3 ) r 2(0 + 1 ) 



g ( r ) = - . (5-12) 

2 2 2 2 2 
(r+y) (r + y ) 

The receiver that implements this test is shown in block 
diagram form in Fig. 34. The ZMNL's for the unit root-mean- 
square Hall models (see Table II) parameter i zed by V D are 
shown in Figs. 35 and 36. The ZMNL receiver char acter l st i c 
for values of V D equal 2.1, 3, 6, 10 and 14 are plotted. 

Recalling that the nonlinearity is optimized for van- 
ishingly small signals the ZMNL character i st l c may be 
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Figure 34. Locally Op-timum Receiver For Energy Detection 



intuitively analyzed. Small values of the sampled envelope 
will be given a negative weighting, intermediate values are 
emphasized and large values of the sampled envelope are 
effectively nulled. As the noise becomes more impulsive, 
which corresponds to a larger value of and smaller 0., the 

above Features become more pronounced . This is seen in Fig. 
35 where the ZMNL has large negative values that extend oFF 
the plot For very small envelope values. The maximum value 
oF the ZMNL occurs at 

1 /2 

x = y( (0+5) / (0+1 ) ) (5-13) 

and the zero crossing at 

1 / 2 

x = y <2/ <0 + 1 ) ) . (5-1 4 ) 

Both oF the points are proportional to the mean value oF the 
underlying distribution For the same Hall model number. 
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Figure 35. Zero Memory Nonl inearity (ZMNL) For Locally Optimu 
Receiver in Hall Noise For = 10, 14 
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Figure 36 . Zero Memory Nonlinearity (ZMNL) For Locally Optimum 
Receiver in Hall Noise For V D =2.1, 3, 6 



Fig. 37 is the P P A versus threshold for the locally 
optimum receiver shown in Fig. 34. The number N of samples 
is equal to 10. The Hall noise samples are par ameter i z ed 
with the same Vo values used in Fig. 32. The threshold data 
obtained from this plot was used to generate the P 0 curves 
shown in Fig. 38. The only complete data set was obtained 
for the Hal 12 case with Vo equal to 2.1. Some data was 
obtained for V 0 equal to 5. The peaks in P 0 for the more 
impulsive noise occur at smaller SNR levels, not shown on 
the plot. The data points are plotted on the graph since 
the dashed curves only approximate the true values. 

It can be seen that the locally optimum receiver does 
improve the P 0 (for the same P F a ) in the small signal case. 
In comparing the values with Fig. 33 for V 0 equal to 2.1 the 
locally optimum receiver has' a performance better than the 
square and sum receiver from -10 dB to 7 dB SNR. Above 7 dB 
SNR, the P 0 for the locally optimum receiver drops to 0. 
This is due to the fact that the receiver is optimized for 
small signals. The Hall2 ZMNL in Fig. 36 shows the reason 
for this. It peaks at ' approx l mate 1 y 1.3 times the RMS value 
of the noise and approaches 0 as the input goes to infinity. 
This behavior of the ZMNL must be modified to obtain a 
receiver that will also work in the large signal case. In 
the next Section a robust receiver structure is proposed 
that does this. A robust receiver is d i f f erent l ated from a 
parametric receiver in that the robust receiver does not 
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Figure 37. P F „ -For Locally Optimum Receiver in Hall Noise 
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depend on the parameters oF the noise density and it per- 



forms well, lF subopt 1 ma 1 ly , over a class or even differing 
classes oF noise densities. 

E. ADAPTIVE LIMITER 

The ZMNL derived in the previous Section is a relatively 
complex function. Furthermore, it does not work For large 
signal levels. In this section we consider a practical 
implementation oF the ZMNL consisting oF an adaptive limiter 
to approximate the ZMNL Function with a modiFication to 
improve large signal perFormance. Although the locally 
optimum receiver was derived For Hall type impulsive noise, 
iF the appr ox l mat l on to it is properly designed, it should 
show robust perFormance in the presence oF^ many impulsive 
noise process. 

Examining the P & curves oF the locally optimum receiver 
For the Hal 12 model shows that it provides increased 
detection perFormance over the square and sum receiver From 
-10 dB to 7 dB SNR. Since this region is where the square 
and sum receiver suFFered the worst perFormance degradation 
in impulsive noise, the Hall2 ZMNL was chosen to be modiFied 
to obtain a robust non 1 i near i ty . The solid curve in Fig. 39 
is the Hal 12 ZMNL. To correct its perFormance at large 
signal levels, the ZMNL was First altered, as shown by the 
dashed line in Fig. 39, by extending its char acter l st i c 
straight out From its maximum value. Furthermore, an 
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equivalent test statistic is obtained if g'(x) is any linear 
function of g(x) [Ref. 46] so 

g'(x) = ocgCx) +6 a ^ 0. (5-15) 




Figure 39. Modification to ZMNL 

With suitable parameters a and 6, the modified ZMNL is shown 
in Fig. 40. The level, a, remains to be determined. 

Recalling Eqn. 5-13, the argument of the ZMNL char- 
acteristic at its maximum value is proportional to the mean 
value of the Hall distribution for which it was derived (see 
Table I). For the Hall2 distribution the proport i ona 1 l ty 
constant is approx l mate 1 y 1.6. Therefore, in order to 

preserve the improvement in P 0 in the same SNR region as the 
Hal 12 ZMNL, the limiter level was set to equal 1.6 times the 
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Figure 40. Linearized Modification to ZMNL 



mean value of the assumed noise distribution, 

A block diagram of the receiver is shown in Fig- 41, 
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Figure 41. Receiver Structure for Adaptive Limiter 
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Note that the ZMNL must be adjusted. This information could 
be known a priori or, as symbolically shown in Fig. 41, 
determined -From the statistics of the noise. Both cases 
will be examined below. 

The false alarm rate performance of this modified 
receiver in Hall and Gaussian noise, where the statistics 
are known a priori, is shown in Fig. 42 for N equal to 10. 
The value for a was determined from the noise mean. At a 
P F a equal to 1 x 10" a , the P 0 curves of this receiver are 
shown in Fig. 43. The threshold was determined from Fig. 
42. The Gaussian noise performance of this receiver is 
poorer than that of the square and sum receiver but only by 
a few percent. However, there is a dramatic improvement in 
the P D for the impulsive noise cases. It is interesting to 
note that for the same noise power, a well designed 
nonlinear receiver provides better detection performance in 
impulsive noise than the optimum receiver for Gaussian noise 
CRef. 49]. 

By plotting the empirically determined threshold from 
Fig. 42 versus the mean of the noise distribution it is seen 
that the relationship is almost linear. This is shown in 
Fig. 44 for two false alarm rates . Using this relationship 
the ZMNL and threshold can be adaptively modified to provide 
nearly constant false alarm rate (CFAR) performance under 
varying noise conditions. This is a desirable charac- 
teristic for an energy detection receiver. 
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Figure 42. P F * For Limiter Receiver in Gauss and Hall Noise 
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Figure 44. Threshold versus Mean Value of Hall Distributions For 

. Limiter Receiver and P F * = 1 x 10- ‘ and 1 x 10- a 



A -Pinal question -to be resolved in implementinq the 
receiver is the real time det erm 1 nat 1 on o-F an estimate of 
the mean and mean square values o-F the noise. One approach 
is to sample the noise in a signal -Free adjacent channel and 
another is to -Form the sample mean o-F the noise using past 
sample values with no signal present [Ref. 50]. We will 
implement the latter approach and test the adaptive receiver 
on actual - signal and noise samples. In sequencing the 
samples two general approaches may be used. The samples may 
be processed in blocks with no overlap which is the case 
correspond 1 ng to the hypothesis tests discussed in this 
Chapter. The other option is sequentially processing each 
sample with rank two updates o-F the test statistic, the ZMNL 
and the threshold. In this case the hypothesis test 
-Formulation will provide a lower bound on P & . 

The sequential sampling approach was chosen and two runs 
o-F the adaptive receiver are shown in Figs. 45 and 46. The 
test statistic, shown by the dotted line, was -Formed -From 
the 10 previous samples using the receiver structure in Fig. 
41. The adaptive threshold, shown by the dashed line, was 
-Formed From the sample mean and mean square estimates oF the 
previous 20 though 50 samples. The sample mean is shown by 
the heavy solid line and the individual envelope data points 
by the light solid line. 

This receiver structure requires H 0 to be true at the 
start oF the test and will detect the transition to H, . 
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Figure 45. Adaptive Limiter Operation in Atmospheric Noise 
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Figure 46. Adaptive Limiter Operation in Gap Noise 



Since -the adaptive threshold assumes noise only samples, it 
will be corrupted when signal plus noise samples are used to 
-Form it. This will occur 20 samples after H t becomes true. 

This effec t can be seen in Fig. 45 where an actual 
signal From one Filter oF a compressive receiver is plotted. 
The noise is highly impulsive and would be a very diFFicult 
detection environment For a receiver optimized For Gaussian 
noise. Samples 50 to 110 are noise only and samples 110 to 
220 are signal plus noise. The test statistic exceeds the 
threshold at sample 125 and the receiver would declare Hi 
true at this point. However, due to the Fact that the test 
statistic and the threshold are both corrupted by the signal 
plus noise samples, the test statistic later drops below the 
new threshold. 

In Fig. 46, the receiver perFormance in a simulated gap 
noise environment is analyzed. The signal commenced at sam- 
ple 110 and the receiver declared the signal present at 
sample 117. Again the threshold is later corrupted and the 
test statistic drops below the threshold. 

F. SUMMARY 

This Chapter has presented a systematic development oF 
the bandpass energy detection problem in impulsive noise. 
Starting with a Hall density Function For the noise, the 
locally optimum receiver is derived. This receiver is 
modiFied to design a practical receiver with robust 
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performance in non -Gauss 1 an noise. 



The design methodology 



clearly delineates the role of the test statistic, threshold 
and nonl 1 near ity . This allows independent evaluation of 
each of these features. One aspect of energy detection 
receivers not considered is the role of data in filters 
adjacent to the filter being processed. In an impulsive 
noise environment, the noise statistics of adjacent filters 
will be correlated and their use may potentially offer a 
further performance gain. 
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VI. CONCLUSION 



A. RESULTS 

The dissertation has investigated the modeling o-P 
atmospheric and man-made radio noise and applied -the results 
-to -the bandpass energy detection problem. 

A generalized means o-P describing an impulsive noise 
process in -terms o-P i-ts complex envelope was derived. 
Previous models o-P atmospheric noise were surveyed and some 
extensions -to -the Hall model for atmospheric noise were 
deve 1 oped . 

Using -the physical characteristics o-P gap noise, a nine 
parameter model was constructed -that: allows arbitrary sour- 
ces o-P gap noise -to be synthesized. The driving impulses 
for -the gap noise were postulated -to be points o-P a 
branching renewal process. The application o-P -this "type o-P 
probabilistic model -to power line noise is unique. 
Additional results include a new proof o-P the power spectrum 
o-P an equilibrium renewal point process and a derivation o-P 
the power spectrum o-P a truncated branching renewal point 
process. The postulated model was shown to correctly 
predict all the significant -features o-P the observed 
spectrum . 

In a similar -fashion, corona noise was postulated as a 
-filtered impulse process with a periodically modulated rate 
-function. This model was shown to correctly duplicate the 



128 



observed time domain behavior of corona noise. A method for 



estimating the parameters of the model was discussed and an 
actual example o-P corona noise simulated. 

The energy detection problem was also considered and 
starting with an assumed density for impulsive noise, a 
design methodology for a parametric and robust receiver was 
developed. This receiver was then tested against actual and 
simulated data and shown to be superior to the Gaussian 
noise receiver when corrupted with impulsive atmospheric or 
power line noise. 

It has been stated that in engineering design one seeks 
not so much to be optimum but to avoid crippling non- 
optimalities CRef. 51J. It is hoped that this research will 
allow system designers to test their systems with simulated 
interference based on the models presented. The robustness 
of the system in an actual noise environment can then be 
evaluated, as was done in this research for the locally 
optimum energy detector. The noise models developed here 
should be particularly applicable to systems that operate 
from a fixed site within line of sight of power lines where 
the chance of having interference from gap or corona noise 
sources is significant. 

B. FURTHER RESEARCH 

At the conclusion of this dissertation several problems 
are worthy of further study. Research should be done to 
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de-Pine the range o-F values for the power line noise model 
parameters and to determine the e-F-Fect o-F power line 
construction practices on the parameters. The suscept- 
ibility o-F various types o-F communication equipment to 
simulated gap and corona noise should also be examined. The 
corona noise model needs to be extended so that it can be 
specified and estimated using a bandpass vice lowpass f 1 Iter 
impulse response. The work on energy detection suggests 
many opportunities -For -Further research. The problem o-F 
estimating the unknown noise parameters -For the adaptive 
updates and the inclusion o-F adjacent -Filter samples into an 
algorithm are prime examples. 
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APPENDIX A 



INSTRUMENTATION 



The instrumentation configuration employed to provide 
data on the detailed time- and frequency-domain properties 
of 2 high -Frequency (HF) radio noise is shown in Fig. 47. The 



£ FIELD 




Figure 47. SimpliFied Block Diagram o-F Measurement System 

radio Frequency CRF) input to the 1 nstrumentat 1 on was either 
a Fixed HF long wire antenna or a whip antenna mounted on a 
mobile van. One complete measurement system is trans- 
portable and may be moved From site to site. The other 
system is installed in a mobile van For noise measurements 
directly at the interFerence source CReF. 533. 

In the lower halF oF Fig. 47, a Hew 1 ett -Packard 141T 
Spectrum Analyzer is used as a scanning receiver to drive a 
Develco 7200B 3- Axis Display. The spectrum analyzer can be 
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tuned to any desired -Frequency in the HF band. Its scan 
rate, scan width, intermediate -Frequency (IF) bandwidth, IF 
gain, RF attenuation, and other controls can be adjusted to 
best describe the noise under observation. 

An alternate and complementary part of the system is 
used to examine the narrowband properties of HF noise and is 
depicted in the upper hal-F o-F Fig. 47. An HF receiver is 
used as an ampl if ler/translator and tuned to a Frequency 
where the noise is present. The demodulated audio output is 
then applied to a Wavetek (Nicolet) UA500A spectrum analyzer 
which subdivides the audio-output spectrum into about 500 
Frequency segments. When the Full audio-output bandwidth oF 
5 kHz is examined, the UA500A provides a Frequency 
resolution oF 10 Hz per segment. The Wavetek analyzer can 
provide individual transForms or average a selected number 
oF successive transForms. The transForms (either individual 
or averaged) are then presented on the 3-axis display or an 
osc i 1 1 o scope . 

The 3-axis display provides a continuous moving real- 
time visual representat i on oF the analog output From the 
spectrum analyzers or narrowband receiver audio output. The 
receiver or analyzer output analog data are digitized in the 
display and stored in a semiconductor memory. The data in 
the memory are Formatted and shown on a cathode ray tube 
(CRT) in a convenient Frequency-ampl ltude-time (3-axis) 
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-Formal:. The 3-axis presentation can be frozen at: any desired 
time and photographed with a standard oscilloscope camera. 

Fig. 48 illustrates the procedure used to make the 3- 
axis present at 1 on . The analog input is divided into 512 




FREQUENCY 

and/or 

SCAN TIME AXIS 



Figure 48. Diagram of Data Format 

equally spaced data points (indicated by the horizontal dots 
of Fig. 48). The signal amplitude at each data point is 
represented by an 8-bit word. When a scan is completed, its 
data are also stored in memory. Line 1 in the view moves to 
line 2, and the new scan appears as line 1. Subsequent 
scans move earlier data, line by line, upwards along the 
time axis to create a rising raster type display. When the 
memory is full (64 scans), each new scan is entered into the 
bottom line, and the oldest data at the top line is 
discarded. The resulting animated view provides a visual 
picture of noise and signals within the block of frequencies 
if the spectrum analyzer output is being observed. If the 
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receiver audio output is being observed the display is a 
stacked series of consecutive time records, 

The 3-axis display has a number of controls to aid the 
observer in analyzing the structure of data presented. 
Among these controls are: Cl) elevation and azimuth geometry 
controls to vary the viewing aspect, (2) an amplitude 
compression control, (3) a threshold control to vary 
background levels, (4) time-axis expansion controls, and (5) 
a stop-action switch to -Freeze the data in memory for 
detailed observation, for observation from various aspects, 
or for photographing. These controls are used to optimize 
the presentation of desired signal detail. 

Accurate frequency , time and amplitude calibrations are 
maintained so that the resulting 3-axis photos can be 
manually scaled for precise signal detail. The digitized 
data at the display input are available at a digital 
interface connector for external digital recording or 
processing. Received data can be processed by computer to 
obtain conventional noise statistics. This feature permits 
the comparison of selected 3-axis views of noise with 
standard statistical descriptors of noise. 
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APPENDIX B 



FILTERED IMPULSE PROCESS 

The complex envelope of a noise process can be described 
by -the Following equation 



Net) -jet 

n(t) = gCt) + E a x e h(t-t t ) 

i = l 


(B- 1 ) 


where g(t) is a quadrature Gaussian process, hCt) 


is a 



linear f 1 Iter, a* is an independent:, identically distributed 
random amplitude and 0 t is a random phase uniformly 
distributed over 0 to 2jt . The joint character i st i c Function 
oF the inphase and quadrature components oF the above 
process is 



j (ait n c (t) + ti) a n a C t ) ) 

0(u)i,u) a )=ECe J. 


(B-2) 



IF the Gaussian noise is assumed independent oF the 
impulsive noise the character i st 1 c Function oF the Gaussian 
term alone can be written 



j (u)i g c Ct) + u>a g a (t) ) 

0 Q ( u)i , u) a ) = E C e 1 . 


( B - 3 ) 


The inphase and quadrature components oF bandpass 
Gaussian noise are independent so 


wh i te 


2 2 

-No B ( U)i + Q)a ) /2 

0 S C u)i , u) a ) = e 


(B-4) 


The joint character i st i c Function oF the ith pulse 

inphase and quadrature impulsive noise terms is 


oF the 
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ja, <(Ji (h c (t-t 4 )cos0i - h fl (t-t 4 )sin0 4 ) + 

0 1 ( tu i , ojg ) = ECe 

(j 2 (he ( t-t 4 )sin0, + h 3 (t-t t )cos0, )) 

: . (B-5) 

Since "the noise process is narrowband, its joint density 
-Function is circularly symmetric and the -Following 
transformation can be applied 



2 1/2 
(Ja ) 



-1 



(Jr = ((Jt + (Ja ) a = tan ((j t /(Jg ) 

(Ji = (j r cosa (j a = (JrSina. (B-6) 

This type o-F transformation o-F a 2-D Fourier transform to a 
one dimensional transform of the envelope is known as a 
Hankel Transform and the inversion -Formula is 



p ( r ) = 



r (Jr 0 ((Jr ) Jo ( r(Jr ) d(Jr 



(B-7) 



0 



where J 0 is the ordinary Bessel Function of the First kind, 
order 0. To simpl lFy the derivation it is assumed that 
h 9 (t) is zero (see comment beFore Eqn . 2-21 ) . Now 

transForming Eqn. B-5 and taking the expected value oF 0 4 
and a ( 

co 2 n 

' 1 j(j r a, h c (t-ti )cos(0i -a) 

0 1 ( (Jr ) = E C p(0t)p(a 4 )e d0 t da, ] 

j t 

0 0 (B-8) 

Since 0, is uniformly distributed over 0 to 2n , 
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OD 



01 ( CD. ) 



EC 

t 

0 



p ( a t ) Jo ( cd. a i h c (t ) ) da i 1 . 



(B-9) 



I -P N(T) in Eqn. B-i is a Poisson process with rate X and 
observation interval T, the T t 's are independent and uniform 
over (0, T) allowing the character l st l c function of the 
impulsive term to be written as 

m k 

0, (cd. ) = E 0 (co r ) PrCN = kD. (B-10) 

k = 0 



Now using the probability law for k events of a Poisson 
process 

a> k -XT k 

0! (ID. ) = E 0 (CD.) e (XT) /k! (B-ll) 

k=0 



jhich sums to 



p ( a ) 



CJ 0 (cd. ah c (t) ) - 1 ] dT da 



(B- 12) 



0 0 

0i (cd. ) = e 

Using the convolution property of zero order Hankel 
transforms CRef. 54] and Eqn. B-7 the density function of 
the envelope for Gaussian plus impulsive noise is 



p ( r ) 



CD 

riD, Jo ( ru),. )0 a ((Ur )0i (w,. ) dw,. . 



~a> 



(B-13) 
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APPENDIX C 



SIGNAL PLUS NOISE DENSITY 

For the truncated Hall model with 0 = 3, the joint 
density Function of the in-phase and quadrature components 
is given by: 



2 2 



p(x,y) 



D 2y 



2 2 2 2 2 



2 2 2 

0 < x + y < Tp 



where 



2n(D - 1 ) ( x + y + y ) 
0 



2 2 2 

x + y > T P 



2 2 1/2 



CC- 1 ) 



D=(l+T P /y) 

The sum of a random phase sinusoidal signal of amplitude A 
and noise can be expressed in terms of its in-phase and 
quadrature components w and z where: 

w = Acos(i|i) + x 

z = Asin(\|/) + y. (C-2) 



Eqn. C-2 is substituted into Eqn. C-l and transformed to 
polar coordinates with 

w = rcos(0) 

z = rs l n ( 0 ) (C-3) 

where r represents the envelope and 0 the phase of the 
signal plus noise. The joint probability density Function 
p(r,0) conditioned on \\j is 
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2 



2 



D r2y 

p(r, 0 |\y) = ... 

2 2 2 2 2 

2n(D -1 )(r - 2Arcos(i \j - 0) + A + y ) 

2 2 2 
0 < r - 2 Arco 5 (i|/ - 0) + A < T P 

= 0 elsewhere. (C-4) 

This expression is valid -For all values o-F r between 0 and 
T P - A. Since our principal interest is the small signal 
case, A << Tp , this will not a-F-Fect subsequent analysis. By 
integrating Eqn. C-4 with respect to 0, the conditional mar- 
ginal density -Function o-F the envelope, p(r| 4 /), is obtained; 

n 



p C r | 4 /) 



p ( r , 0 | 4 /) d 0 



0 



CC-5) 



The integral is a -Function only o-F cos ( 4 / - 0 ) , 4 / 

nowhere else in the integral. Therefore by setting 
0 and using Eqn. 2.554-3 in CRe-F. 543 gives 

2 2 2 2 2 



p C r | 4 /) 



( r + A + y ) 2y rD 



4 22 22 22442 

( r - 2A r + 2A y + 2y r + A + y ) (D - 1) 



appears 
a = 4 / - 



(C- 6 ) 



Since 4 / is assumed independent -From sample to sample and 
uniformly distributed over 2 n , 

2 n 



p (r ) 



p ( r | 4 /) p ( 4 /) d 4 / 



( C - 7 ) 



and 



0 
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2 2 2 2 2 
( r + A + y ) 2y rD 



p ( r ) = 



4 22 22 22 442 

(r - 2A r + 2 A y + 2y r 2+ A + y )(D - 1) 











(08) 


By a similar argument, 


■the density -Function 


-For a 


random 


phase sinusoidal 


s i gna 1 


in narrowband noise 


with a 


Hal 15 


distribution is 


p ( r ) = 


4 2 2 

2y r ( b + 2a ) 




(09) 



2 2 5/2 

(a - b ) 



where 



2 2 
a = r + A 



and 



b = -2Ar. 
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APPENDIX D 



SPECTRUM OF A RENEWAL PROCESS 



The complex envelope representation of 
noise on a Finite interval T is given by 

N ( T ) -j0 4 

n(t) = E a 4 e h(t-t 4 ) ; 

i = 1 



Filtered impulse 



0 < t < T. 
(D- 1 ) 



N(T) is a random unit counting process denoting the number 
oF impulses in the interval, the -Ct 4 > are the random arrival 
times and 0 4 is uniFormly distributed over 0 to 2n . The 
pulse amplitudes -Ca 4 > are identically distributed 

stat i st i ca 1 1 y independent random variables with second 
moment A 2 and Fourth moment A 4 . The Filtering eFFect on the 
impulses are accounted For by the time invariant complex 
envelope impulse response h(t). 

The envelope squared oF n(t) is 

2 N (T) 2 

E (t) = n(t)n" (t) = E a 4 h(t-t 4 ) h" (t-t 4 ) 

i = l 



( D-2 ) 

where we have assumed terms oF the Form h(t-t 4 ) h" (t-t k ) For 
di FFerent i and k are negligible. We estimate the power 
spectral density (PSD) oF the envelope by averaging 
magnitude squared , length T, Fourier transForms oF E a (t). 
The mean value oF the estimate is 



S ( oi , T ) = ECS* <iD)e a " <w) D/T 



(D-3) 
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where 



£ a <oj,T) 



T 

2 - jiut N(T) 2 -j(ut 4 

E (t)e dt = E a t e H a (w) 



and 



H a (oj) 



-jwt 

h(t) h* (t)e dt 



u 

CD 

is the Fourier transform of the magnitude squared impulse 
response- We have assumed the impulse response is much 
shorter than the observation interval T. 

The mean o F the estimate becomes 



S(cu,T) 



N(T) N(T) 2 2 -joj(t t -t k ) 2 

EC E E a t a k e ] | H a (w) | /T. 

k=l i=l 

( D - 4 ) 



This is a compact and general expression tor the PSD o-P a 

truncated Filtered point process. It is in terms o-P a 

summation of the charact er i st i c Functions oF the 

interarrival time distributions between all combinations oF 

pairs oF points. (The term i ncor por at i ng the Frequency 

response oF the Fi Iter will be set to one For the rest oF 

the analysis.) Conditioning on N(T) 

_ cd N N 2 2 - j (i) ( 1 1 -t k ) 

S((d,T) = E EC E E a t a k e DPrCN = 1]/T. 

1=2 k= 1 i=l 



CD-5) 

Plotting the diFFerence k-j indicates a diFFerent summation 
over m and n shown below. 
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k N n N 



-3 


-2 


-1 


0 






-3 -2 -1 


1 


2 3 


-2 


-1 


0 


i 






-2 -1 


1 


2 


-1 


0 


1 


2 


. . . 


N 


-1 


1 




0 


1 


2 


3 




J 









Let m = i - k and n = i + k such that t« = t 4 - t k is the 
mth interarrival time. Then Eqn. D-5 can be written as 



_ oo N-l -n 2 - j wt, 

S ( tu , T ) = E E-C £ E A a e + NA 4 

1=2 n= 1 m=l 

N- 1 n 2 . j tut. 

+ £ E A a e >PrCN = 1J/T. 

n= 1 m=l ( D -6 ) 



On -Further rearrangement 

— _ oo 2 N-l n j cut. 

S ( 0 ) , T ) = NA 4 /T + E A a E-C E E e + 

1=2 n=l m= 1 



- jcut. 

e >PrCN = 1 ]/T 

CD-7 ) 



and conditioning the interarrival time distribution over the 
number o-P pulses N in the interval T 



S ( tu , T ) 



where 



_ 2m 1-1 n 

NA 4 /T + Aa E E E [0. (j(u|N=l) + 

1=2 n= 1 m= 1 

0. " ( j tu | N=1 ) ]Pr[N = 13/T 

(D-8) 



j(Ut. 

0 m ( j co | N= 1 ) = ECe 1 

is the character i st i c -Function o-P the m* h interpulse spacing 
given exactly 1 pulses in the interval. The -Following 

assumption is now made; that the mth interpulse spacing is 
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■the sum o-f m statistical ly independent interpulse sep- 
arations, and is independent of 1, the number o-f pulses in 
the interval. This assumption will give an asymptotically 
correct result tor the PSD as the observation interval T, 
becomes much greater than p, the mean value o-f the 

With the above assumption 

m 

0. Cj(i)|N=l) = 0 (JOJ) CD-9) 

where 0(jid) is the char acter i st i c -function o-f the inter- 
pulse separations. Eqn. D-8 then sums to 

n+ 1 

_ _ 2 a) 1-1 0 ( j a) ) -0 ( j a) ) 

S ( a) , T ) = na 4 /t + A a E EC — — — + 

1=2 n= 1 1 -0 C j a) ) 



n+1 



0" ( j a) ) -0" C j a) ) 
1 -0* (JO)) 



UPrCN = 1 ] / T 



CD- 10) 



and summing again 



1 + 1 



„ _ 2 CD 10( j(D) 

S(w,T) = NA 4 /T + A 2 E -C 

1=2 



0 ( j 0) ) 



0 ( j ID ) 



2 2 

(l-0(j (D)) (l-0(jtD)) ( 1 —0 ( j tD ) ) 



1 + 1 



10* ( j(D ) 



0" ( j(D ) 



0" ( j 0) ) 



2 2 

(1-0* ( j (D ) ) C1-0*(J(D)) (l-0*(jCD)) 

Rearr ang l ng 

_ 2 CD 0 ( J (D ) 

S ( (D , T ) = NA 4 /T + A 2 /T<N - PrCN=i:> E C 

1=2 



>PrCN = 1II/T. 



(D-l 1 ) 



0" ( JO)) 



( 1 -0 ( j (D) ) C 1 -0 ( j(D) ) 
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2 0CjUl) 

A a /T <1 - PrCN=0]-PrCN=l> < 



0 * C j id ) 



2 2 
( 1 -0 ( joi) ) (1-0- C j(il) ) 



1 + 1 

2 00 0 ( j li) ) 

+ A a /T E PrCN = 1 ] < 

1=2 2 

( 1-0 ( j(D) ) 

Since PrCN = 1] is less “than 1 
41:p. 115] -For any distribution 

convergence o-P -the -third term is 
equilibrium renewal process 



1 + 1 

0* ( jlD) 

> . 

2 

( 1 - 0 - ( j (D) ) ( D - 1 2 ) 

and |0 ((d)| < 1 C Ref - 

then -the absolute 
guaranteed . For an 



N = ECN(t) ] = T/p CD-12) 
where p is the mean time between renewals CRe-P. 32 : p . 46]. 
Now letting T, the observation interval, go to in-finity the 
asymptotic PSD is 



__ A a A 4 0 ( J ( 1 ) ) 

SC(D) = lim SCui,T) = < + 2 ReC ] > . 

T = >00 p 2 1-0 ( j a) ) 

A a CD-13) 

This is the expression -For the PSD o-P an independent 

increment point process in terms o-P the character 1 st 1 c 

-Function o-P the interarrival distribution. A similar result 

result was developed by Cox and Miller by considering the 

renewal intensity -Function CRe-P. 55]. 
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APPENDIX E 



SPECTRUM OF A BRANCHING RENEWAL POINT PROCESS 

The complex envelope representation of f 1 ltered impulse 
noise arising -Prom a branching renewal process on a finite 
interval T is given by 

M(T) N. - j 0 ■ n 

n(t) = E E a. „ e h (t“t. „ -T. ) CE-1) 

m=l n=l 

where T. is “the beginning of the mth main interval, N. is 

the number o-P impulses in the mth interval, arriving at 

times -Ct. „ > after its onset and M(T) is the unit counting 

process defining the number o-P main events in the 

observation interval T. The pulse amplitudes -Ca. „ > and 

phases -C0.„> are stat i st i ca 1 1 y independent random variables 

and the -C0.„> are taken to be uniform on -C0,2n>. The net 

Filtering effects on the impulses are accounted For by the 

time invariant complex envelope impulse response, hCt). 

The envelope squared oF the process is 

2 M(T) N. 2 2 

E (t) = n(t)n*(t) = E E a . „ | h ( t-t. n + T. )| CE-2) 

m=l n= 1 

and using the same assumptions as Appendix D the mean value 
oF the Bartlett estimate oF the PSD is 

2 

S ( (d , T ) = EC | £* (u>,T) | D/T CE-3) 

where ( co , T ) is the length T Fourier TransForm oF E 2 (t) . 

Spec l Fi ca 1 ly , 
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M(T) N. 2 

£ a ( cd , T ) = E E a * „ e 



- jut. „ - j(dT. 



e 



Ha ( Id ) 



(E-4) 



m=l n= 1 



T 



With Ha (cd) 



2 -jcdt 

|h(t)| e dt , -the Fourier transform of the 



t 



0 



magnitude squared impulse response. To Find the mean oF the 

estimate we must evaluate 

_ 2 M(T) N. M(T) Nj 

S((d,T) = EC | H a (td) | E E E E 

m=l n=l 1=1 k=l 

2 2 ” j (d ^ t B n ”t i k ) “jcdCT. - Ti ) 

a. n a 1 k e e 3 /T 

<E-5) 

In order to simpl lFy this expression, we make use oF 
assumptions 1, 2 and 4 From section 3D; (1) ECa 2 .„a 2 i k 3 = 

a 4 , which will be normalized to one For simplicity, (2) N m 
and Ni equal either N1 or N2 alternating between the two 
values and (3) F(t.„) and F(ti„), the interarrival time 

d i str i but i ons , alternate between the negative and positive 
phase distribution parameters . Furthermore, let M be an 
even number oF main process points occurring at intervals, 
T 0 /2, oF a Fundamental Frequency (in our case F 0 = cd 0 /2n = 

1/To = 60 Hz, the power-line Frequency) and let M2 = M/2. 

This assumption will greatly simpliFy the calculations and 
For M greater than 10, not s i gn i Fi cant ly reduce the accuracy 
oF the result. 
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Now looking 



only at the factors involved in the 



summation in Eqn. E-5, six separate terms can be identified: 



Term 1 - The inter-pulse group summations between pulses in 

pulse groups with N1 pulses 

N 1 N 1 jo) (t„ -t k ) 

= M2 E E ECe 3 

n= 1 k= 1 

( E - 6a ) 

Term 2 - The inter-pulse group summations between pulses in 

pulse groups with N2 pulses 

N2 N2 j (D ( t n -t„ ) 

= M2 E E ECe 3 

n= 1 k= 1 

(E-6b ) 

Term 3 - The intra-pulse group summations between pulses in 

groups with N1 pulses 



M2 


M2 - j (D ( m- 1 ) T 0 


N 1 


N 1 jut n -JW'tk 


E 


E e 


E 


E ECe ] E C e : 


m= 1 


1 = 1 


n= 1 


k=l 



m< > 1 



( E - 6c ) 



Term 4 - The intra-pulse group summations between pulses in 
groups with N2 pulses 





M2 M2 


-jd)(m-l )T C 


N2 


N2 


jwt n 


“ JCDth 




= 


E E 


e 


E 


E 


ECe 


] E C e 3 






m= 1 1=1 




n= 1 


k= 1 










m< > 1 










<E- 


6d ) 


T erm 5 


- The intra-pulse 


group 


summat l ons 


between pulses 


l n 


groups 


with N 1 


pulses per 


group 


to 


pu 1 ses 


in groups with 


N2 



pu 1 ses 
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j(D<To /2-T D ) M2 M2 -j(i)(m-l)To N1 N2 jwt„ -jut* 



= e 


E E e EE ECe lECe 3 

m — 1 1 = 1 n=l k=l 

(E-6e) 



Term 6 - The intra-pulse group summations between pulses in 

groups with N2 pulses per group to pulses in groups with N1 
pu 1 ses 

- ju)(T 0 /2-To ) M2 M2 - j cu ( m- 1 ) T 0 N2 N1 ju)t„ -ju>t k 



= e 


E E e EE ECe lECe ] 

m=l 1=1 n=l k= 1 

(E-6f ) 


It is 


important "to note the expression -For the 



character i st i c -Function o-F the interarrival times differs 



between the 


inter-pulse and intra-pulse cases. In the 


i nter-pu 1 se 


case, discussed in Appendix D, the summation is 


only over 


a -Function o-F the difference in the pulse 



positions where in the intra-pulse case the summation is 
over the absolute position o-F each pulse in its group. 

Terms 1 and 2 can be evaluated using the results -From 
Appendix D and assumption 3 -From section 3D: 



Term 1 = 


M2 5(ta, Nl) (E-7 a) 


Term 2 = 


M2 5(0), N2) <E-7b) 



and the remaining terms can be determined by straightforward 
evaluation o-F the -Finite sums 



Term 3 = 


cos(M2 2no)/o)„ ) - 1 

•C - M2 > fi(o),Nl) 

cos ( 2n 0)/u)o ) - 1 (E-7c) 


Term 4 = 


cos(M2 2ntD/(D 0 ) - 1 

{ : - M2 > (1(0), N2) 

co5(2n(o/(i) 0 ) - 1 (E-7d) 
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Terms 5 and 6 are complex conjugates and can be combined to 
•Form one term 



cos(M2 2ncu/u) 0 ) - 1 

Term 5a = 2cos (a)(T 0 /2-T & ) ) { — — — > ( 0 1 , N1 ,0 a , N2) 

cos (2ntd/td 0 ) - 1 



C E-7e ) 

where 0i and 0 a are the char acter i st i c -Functions associated 
with the interarrival times -For each phase and N 1 and N2 are 
the number o-F renewals associated with each phase. The 
-Functions ft, and i \s are de-Fined below; 



§ <0 , N) 



2 

N - N | 0 ( j a) ) | 



2 

| 1 - 0(jio) I 



2 Re< 



N 

0 ( jCD) ( 1 - 0 ( jU) ) 



2 

( 1 - 0 ( jtu) ) 



>, 



(E-8) 



and 



Q C0 ,N) 



N+l 2 

|0(jlD) - 0(j(l)) | 



2 

I 1 - 0( JO)) I 



(E-9) 



( 0 1 , N 1 ,0 a ,N2) 



Nl+1 N2+ 1 

(0 4 (jcd) “ 0i (jo)) ) (0 a (-jo)) - 0 a ( - j a) ) ) 



( 1 - 0, (JO)) ) (1 - 0 a ( - j(D) ) 

( E - 1 0 ) 



Thus 



5 

S ( id , T ) = -CECH a ((d) 1/jy E Term(n). 

n= 1 



(E-l 1 ) 
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APPENDIX F 



NON-HOMOGENEOUS POISSON PROCESS 

The -Following derivation is -For the mean and variance of 
a compound periodic non- homogeneous -Filtered Poisson point 
process. This proof is a gener a 1 1 z at 1 on o-F a proof 

presented in Papoulis CRe-F. 4 1 : pp . 382-383D. The expression 
-For such a process on an interval 0 to T is given by 

N ( T ) 

n(t) = E a t h(t-ti ) , (F-l) 

i = l 

where NCT) is a unit counting process with periodic rate 
parameter X(t) , a t • is the weighting o-F the ith point and is 
independent and identically distributed (I ID) -From point to 
point and hCt) is the time invariant impulse response o-F the 
Pi Iter. 

The time axis is divided into consecutive intervals Ij 
of length aoc where aoc = ocj - 0 Cj - t . The number of jumps in 
the counting process in the j t h interval Ij is given by Am a . 
IP aoc is suPPiciently small, then the contribution to the 
total noise process Prom the j t h interval is 

Am 4 

Anj * E a t h(t- 0 Cj ) , (F-2) 

1 = 1 

where Amj is a Poisson random variable with its rate 
parameter approx i mate 1 y equal to 

XCoCj ) a oc • CF-3) 
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The moment function oF the 



interval is then 



j 



h 



Arrij 

s E a* h ( t -oCj ) 
i = l 

a0 j (s) * ECe : . ( F -4 ) 

Conditioning on Arrij , 

ArOj 

s E a t h(t-cXj) 
i = l 

a0j Cs) * ECe | Airij = k] PrCAirij = k] CF-5) 



and 



k sa t h ( t “OCj ) 

a0j Cs) * EC n e ] PrCAmj = kD. CF-6) 

1 = 1 

Since the a t 's are independent -Prom pulse to pulse 

k sah ( t -oCj ) 

as j Cs) * n ECe ] PrCArHj = k], CF-7) 

1 = 1 



where the expectation is with respect to the random 
amplitude a. Using the Poisson probability law and recog- 
nizing the series expansion for the exponential gives 



sah C t-a j ) 

XCttj )Aa-CECe ] 

A0 j C s ) * e 



S i nee 

ao 

n C t ) = E An, C t ) , 

1 = “CD 

a sum oF independent random variables, then 



1 > 



CF-8) 



CF-9) 



GO 

0 n Cs) = n A01 Cs) CF-10) 

1 = -0D 

using the convolution property o-P the moment generating 
Function. The cumulant generating Function is deFined as 
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H/„ (s) 



lnC0n (s) 3 



(F-l 1 ) 



where In 
As Aa => 



Using a 
funct i on 

V|/ n (Si , S 8 



is *the natural logarithm and 



H/ n (s) = E 1 nC A0 t ( s) 3 . 

i = -OD 



0, then 



( s ) = 



p ( a ) 



sah ( t-oc ) 

X ( a ) C e -1 3 dada 



(F-12) 



( F - 1 3 ) 



similar argument the joint cumulant generating 
•For two random variables, n(t t ) and n(t a ) is 



CD 



CD 






p ( a ) 



J L 

~CD -CD 



X ( a ) 



Si ah ( t t -a ) 

C e 



s a ah ( t a -a ) 

- 13 dada . 

( F- 1 4 ) 
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APPENDIX G 



AUTOCORRELATION OF CORONA NOISE 



The mean value o-F -the cyclostationary process defined by 
assumptions 1 and 2 in section 4D is given by Eqn. 4-12. 
Substituting Eqn. 4-20 gives 

OD 2 



ECe(t) ] 



ECa] 



N t 1/2 -as 

E a„ cos (<u# n (t-s) ) (a/n ) e ds . (G-l) 

n=0 



0 



Expanding the cosine and integrating 

2 2 

• - (Uo n 



Nt 

ECe(t)] = (Bn/ 4) E C e 

n=0 



4a 






(Uo n a> 

E 

1/2 k = 1 

2 ( an ) 



1 



( 2k- 1 ) « ! 



1 /2 



a„ cos ( (Uo nt ) + 



2 2 

-(Uo n k-1 

( ) a„ sin((Uo nt) ] 

2a 

(G-2) 



where (Bn/2) = ECa] and (2k-l) 1 ! =135.. .2k-l [Ret. 54: 

Eqn. 3.897]. The sine term is in quadrature with the 

dominant cosine term o-F Eqn. G-2 and represents a phase 
shi-Ft in the mean o-F the process relative to the rate 
-Function. This is due to the -Filtering effect where the 
impulse response o-F the -Fi Iter persists, thus causing a lag 
in the mean value o-F the process. For the values o-F a and 
(Uo considered in this work, the sine term in Eqn. G-2 will 
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be negligible relative to the cosine term. Physically, this 
means that -For -Filter bandwidths much greater than the 
-Fundamental -Frequency o-F the corona noise, the phase shift 
will be insignificant. 

Similarly the variance o-F the envelope when considered 
as a cyclostationary process is 

2 2 
- ido n 



VarCe(t)] * 213 



N t 1 1/2 8a 

E C — ( cc /2n ) e 

n=0 2 



a n cos ( o)o nt ) D . 

(G-3) 



Using Eqn. 4-14 the covariance o-F the assumed process is 



Cov ( ti , t a ) = 213 



cd 2 2 

1 N t -a (ti - s ) -a(t 2 -s) 

E a n cos((i) 0 ns) (a/n )e e ds 

n=0 

max ( t 4 , t 2 ) ( G-4 ) 



where after rearranging and completing the square 



-a (t t -t a ) /2 

Cov(t t ,t s ) = (2f3(a/n)e ) x 



'• N t -2a (s- (ti + 1 a ) /2 ) 

E a n cos((i) 0 ns)e ds . 

J n=0 

max ( ti , t a ) (G-5) 



Lett i ng 

u = (s - (t* + ta ) / 2 ) , (G-6) 

the covariance is 



155 



Cov C“ti , t 2 ) 



-a (-tt 

= ( ( a / n ) 28 e 



2 

-t 2 ) / 2 

) 



x 



CD 2 

1 Nt -2a u 

E a„cos((i) 0 nCu+Ct! +t 2 )/2)e du 
v n=0 

| t t -t 2 | ( G - 7 ) 

which cannot be evaluated in closed form. 

The stationary autocorre 1 at i on function is -Pound by 
using Eqns . G-2 and G-7 and the definitions given by Eqns. 
4-15 and 4-19. The term due to the covariance in the auto- 
correlation -Function is 



2 

-an- /2 

28(a/n)e /T 0 



Nj 

E C a n cos ( Uo n ( 2t ) ) 
n = 0 



-2au 

cos ( o)o nu ) e du 



0 



I ** I 



+ a n sin(a)on(t+ r T/2)) 



2 

-2au 

sin((u 0 nu)e duDdt. 



The sinusoidal terms will integrate to 
exception o-P the a Q cosine term and the 
1 l -P l e s to 

2 oo 2 



-a*r /2 

215 (a /n ) e a 0 



-2au 

e du 



(G-3) 

zero with the 
expression simp- 



(G-9) 



This term will be dominated by the f irst exponential and is 
approximated by 
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2 

-OCT /2 



1/2 



2Be a 0 ( a /8n ) 

The stationary autocorrelation function 
cyclostationary mean value -Function is 



1 



To 



To 



0 



ECe(t) ]ECe(t+T ) D die, 



(G-10) 
term due "to 



(G-l 1 ) 



which using Eqn. G-2 is equal "to 

2 2 
- Uo n 



Nt 4a 2 2 

26 (1/8) E e a n cosoj 0 ni + 2fi(l/4)a 0 . 

n = 1 

CG- 12) 



An approximate expression for the stationary autocorre 1 at l on 
■Function is 



2 

-aT /2 1/2 2 

R(t) % 26Ce a 0 (a/8n) + (l/4)a 0 



2 2 
-ojo n 



N t 

+ (1/8) E e 
n= 1 



4a 2 

a„ cos ( uj 0 nT ) 3 . 

(G- 13) 
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APPENDIX H 



LOCALLY OPTIMUM ENERGY RECEIVER 



The observations consist o-f a sequence o-P N complex 
samples. The hypotheses are: 

H 0 : x t = n 4 

versus 

1/2 

Hj : x t = n 4 + A s 4 i = 1 , . . . , N ( H - 1 ) 

where n 4 are an independent, identically distributed (I ID) 
complex noise samples with a joint inphase and quadrature 
density p ( n c , n 8 ) . The sequence s t is a complex, zero mean 
signal sequence with a known variance, a t , and A is a real, 
positive number. The 1/2 power o-f A is chosen to -facilitate 
the derivation. 

The generalized likelihood ratio CRe-P. 6:p. 585,2111 to 
test the hypotheses described by Eqn. H-l is 



N 

A(x,A) = n 
i = l 



1/2 1/2 
P ( Xc 1 ~ A Sc i , Xg i -A Sfi 1 ) 

p(x c I , X 8 i ) 



P ( Sc i , Ss t ) ds c l ds a 



o 

"CD ”CD 



(H-2) 



Applying Eqn. 5-8, and d i fferent i at i ng with respect to A, 



N 

t L o < x , A) = E 

i = l 



1/2 1/2 

5c 1 ■op < Xc i - A Sc i , x» , -A 5> I ) / 2 



1/2 



■oxc i p (Xc I , Xs » ) A 
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1/2 1/2 

»s i ’opCxc i “A Sc i >Xg i *A Sg i ) / 2 



1 p ( S c 1 9 Sa 1 ) ds c 1 dSa 1 



1/2 



*0X3 i p(x c i ,X S 1 ) A 



A=>0 „ 
(H-3) 



Le-t-ting A go -to 0 and simultaneously -taking -the expectation 
o*F s c 1 and s 9 t requires L ' hosp i*ta 1 7 s rule -to evaluate “the 
expression [Re-F. 56D . A-Fter -taking derivatives -the test 
s-ta-tis-tic is 



N 

t LO (x) = E 
i = l 



2 2 1/2 1/2 
Sc 1 -o p ( x c 1 “A Sc i j X 9 t - A s s 1 ) / 2 



■OXc 1 pCx c 1 ,x 9 t ) 



2 1/2 1/2 
Sc t Sfi 1 *D p (Xc 1 “A Sc 1 9X3 1 “A S3 1 ) 



T 3 X C i -OXa t P (X C 1 , X 9 1 ) 



2 2 1/2 1/2 
Ss 1 -D p ( Xc t - A Sci ,X 8l -A Sgi )/2 

1 p ( s c 1 9 Sg i ) dsc 1 ds s i 



-OXa t p(Xci ,X 9l ) 



A = > 0 . 
(H-4) 



In general, -the inphase and quadrature noise components will 
be circularly symmetric even -though not necessarily 
s-ta-t 1 s-t i ca 1 ly independent: CRe*F. 41 :p. 133D such -that: 

2 2 1/2 



w i*t h 



•F ( r ) = p((x c + x s ) 
2 2 1/2 



( H-5 ) 



r = ( x c + x 8 



(H-6) 
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Now A may be set to 0 and 



N 

t L o (r ) = E 
i = l 



CD CD 

2 2 2 2 
Set Cri-F // (ri)x C t - r 4 -F 7 ( r t ) x c t +■ r 4 1/2 

C — — — — — - - — — — + 

3 

r 4 -PC r t ) 



s c t Ss t C r t t ' 7 ( r t ) x c t x 9 t - t ' ( r 4 ) x c t x 9 t 3 



r 4 -F( r i ) 



2 2 2 2 
Sst Cr t f ' 7 (r t )x fl j - r 4 -f 7 ( r 4 ) x 9 4 + r 4 ]/2 

— — ]p(s c t ? s 9 t ) ds c 4 ds 9 t . 

3 

r 4 t(r 4 ) ( H - 7 ) 

It Set and s 9 t are assumed to have equal variances and a 
zero covariance the test statistic is now written as 



where 



N 2 

t L o ( r ) = E a i g 4 ( r 4 ) 
i = l 



t 7 7 (r) t 7 (r) 

g ( r ) = — + 



( H-8 ) 



( H-9 ) 



•f(r) rf(r) 

I -f the unknown signal components are assumed to have equal 
variances from sample to sample then an equivalent test 
statistic is 



N 

t 7 l o (r) = E g(r 4 ) 
i = 1 



(H-10) 
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